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MEAN FIELD LIMITS OF THE GROSS-PITAEVSKII AND 
PARABOLIC GINZBURG-LANDAU EQUATIONS 

SYLVIA SERFATY 


Abstract. We prove that in a certain asymptotic regime, solutions of the 
Gross-Pitaevskii equation converge to solutions of the incompressible Euler 
equation, and solutions to the parabolic Ginzburg-Landau equations con¬ 
verge to solutions of a limiting equation which we identify. 

We work in the setting of the whole plane (and possibly the whole three- 
dimensional space in the Gross-Pitaevskii case), in the asymptotic limit 
where e, the characteristic lengthscale of the vortices, tends to 0, and in 
a situation where the number of vortices Ng blows up as e —>■ 0. The 
requirements are that Ng should blow up faster than |log£| in the Gross- 
Pitaevskii case, and at most like |log£| in the parabolic case. Both results 
assume a well-prepared initial condition and regularity of the limiting initial 
data, and use the regularity of the solution to the limiting equations. 

In the case of the parabolic Ginzburg-Landau equation, the limiting 
mean-field dynamical law that we identify coincides with the one proposed 
by Chapman-Rubinstein-Schatzman and E in the regime Ng <C |log£|, but 
not if Ng grows faster. 

keywords: Ginzburg-Landau, Gross-Pitaevskii, vortices, vortex liquids, mean- 
field limit, hydrodynamic limit, Euler equation. 

MSG: 35Q56,35K55,35Q55,35Q31,35Q35 

1. Introduction 

1.1. Problem and background. We are interested in the Gross-Pitaevskii 
equation 

(1.1) idtu = Au + ^(1 — in 

and the parabolic Ginzburg-Landau equation 

(1.2) dttt = Art-F-^(1 — luP) in 

in the plane, and also the three-dimensional version of the Gross-Pitaevskii 
equation 

(1.3) idfU = Au + ^(1 — in 

all in the asymptotic limit e ^ 0. 

These are famous equations of mathematical physics, which arise in the 
modeling of superfluidity, superconductivity, nonlinear optics, etc. The Gross- 
Pitaevskii equation is an important instance of a nonlinear Schrodinger equa¬ 
tion. These equations also come in a version with gauge, more suitable for the 
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modeling of superconductivity, but whose essential mathematical features are 
similar to these, and which we will discuss briefly below. There is also interest 
in the “mixed flow” case, sometimes called complex Ginzburg-Landau equation 

(1.4) (a + ib)dtu = Au + ^(1 — in 

For further reference on these models, one can see e.g. [Ti rrniXKiigs^ . 

In these equations, the unknown function u is complex-valued, and it can 
exhibit vortices, which are zeroes of u with non-zero topological degree, and 
a core size on the order of e. In the plane, these vortices are points, whereas 
in the three-dimensional space they are lines. We are interested in one of the 
main open problems on Ginzburg-Landau dynamics, which is to understand 
the dynamics of vortices in the regime in which their number blows up as 
e —)■ 0. The only available results until now are due to Kurzke and Spirn |KS2] 
in the case of ()1.2p and Jerrard and Spirn |.TS2j in the case of (jl.ip . The latter 
concern a very dilute limit in which grows slower than a power of log |loge| 
as e —?> 0 (more details are given below). 

The dynamics of vortices in (jl.ip and (II.2h was first studied in the case where 
their number N is bounded as e ^ 0 (hence can be assumed to be independent 
of e. It was proven, either in the setting of the whole plane or that of a bounded 
domain, that, for “well-prepared” initial data, after suitable time rescaling, their 
limiting positions obey the law 

(1.5) ^ = (a + JI6)V,IT(ai,...,aw) 

where J is the rotation by 7r/2 in the plane, and W is in the setting of the plane 
the so-called Kirchhoff-Onsager energy 

(1.6) W (oi,..., aw) = —VI didj log |aj — Oj| 

where the dj’s are the degrees of the vortices and are assumed to be initially in 
{1, —1}. In the setting of a bounded domain and prescribed Dirichlet boundary 
data, W is (11.61) with some additional terms accounting for boundary effects. It 
was introduced and derived in that context in the book of Bethuel-Brezis-Helein 
[BBH] where it was called the “renormalized energy”. 

In other words, the vortices move according to the corresponding flow (gra¬ 
dient, Hamiltonian, or mixed) of their limiting interaction energy W. After 
some formal results based on matched asymptotics by Pismen-Rubinstein and 
E in |PR[ lElj . these results were proven in the setting of a bounded domain 
bv Lin jLilj and Jerrard-Soner m in the parabolic case, Colliander-Jerrard 
|GJ11 IGJ2] and Lin-Xin |LX2j with later improvements by Jerrard-Spirn |JSpl| 
in the Schrodinger case, and Kurzke-Melcher-Moser-Spirn |KMMSj in the mixed 
flow case. In the setting of the whole plane, the analogous results were ob¬ 
tained by Lin-Xin [LX1| in the parabolic case, Bethuel-Jerrard-Smets |BJS] in 
the Schrodinger case, and Miot |Mi] in the mixed flow case. A proof based on 
the idea of relating gradient flows and P-convergence was also given in [SS4j . it 
was the initial motivation for the abstract scheme of “P-convergence of gradient 
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flows” introduced there. Generalizations to the case with gauge, pinning terms 
and applied electric field terms were also studied Spl , Sp2t IKSli ITIl IST2] . 

All these results hold for well-prepared data and for as long as the points 
evolving under the dynamical law (II.5p do not collide. In the parabolic case, 
Bethuel-Orlandi-Smets showed in the series of papers [B0S1[[B0S2LIBDS3] how 
to lift the well-prepared condition and handle the difficult issue of collisions and 
extend the dynamical law (II.5p beyond them. Results of a similar nature were 
also obtained in |Selj . 

The expected limiting dynamics of three-dimensional vortex lines under (|1.3I) 
is the binormal flow of a curve, but in contrast to the two-dimensional case there 
are only partial results towards establishing this rigorously 


When the number of points blows up as e —>■ 0, then it is expected that 
the limiting system of ODEs ()1.5p should be replaced by a mean-held evolution 
for the density of vortices, or vorticity. More precisely, for a family of functions 
Us, one introduces the supercurrent js and the vorticity (or Jacobian) /ig of the 
map Us which are dehned via 


(1.7) js := {iUs,Vus) /Tg := curljg, 

where {x, y) stands for the scalar product in C as identihed with via {x, y) = 
^^{xy+yx). The vorticity //g plays the same role as the vorticity in classical huids, 
the only difference being that it is essentially quantized at the e level, as can 
be seen from the asymptotic estimate y,s — as e ^ 0, with {oj} 

the vortices of Ug and di G Z their degrees (these are the so-called Jacobian 
estimates, which we will recall in the course of the paper). 

The mean-field evolution for fi = limg^ocan be guessed to be the 
mean-field limit of (II . as N —)• oo. Proving this essentially amounts to showing 
that the limits e —)• 0 and N ^ oo can be interchanged, which is a delicate 
question. 

In the case of the Gross-Pitaevskii equation dnii-doi), it is well-known that 
the Madelung transform formally yields that the limiting evolution equation 
should be the incompressible Euler equation (for this and related questions, see 
for instance the survey [CDS]). In the case of the parabolic Ginzburg-Landau 
equation, it was proposed, based on heuristic considerations by Chapman- 
Rubinstein-Schatzman [CR.Sj and E |E2j . that the limiting equation should 
be 


(1.8) dty, — div {fiVh) = 0 h = —A 

where y is the limit of the vortex density, assumed to be nonnegative. In fact, 
both papers really derived the equation for possibly signed densities, [CRSj did 
it for the very similar model with gauge in a bounded domain, in which case 
the coupling h = —A~^y is replaced hy h = (—A -|- I)~^y (see also Section 
ll.3.4l belowi. and |E2j treated both situations with and without gauge, also for 
signed densities, without discussing the domain boundary. 

After this model was proposed, existence, weak notions and properties of 
solutions to (II.8p were studied in [LZll IDZl ISV] (see also |MZ] for some related 
models). They depend greatly on the regularity of the initial data y. For y 
a general probability measure, the product yVh does not make sense, and a 
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weak formulation a la Delort [Dej must be used; also uniqueness of solutions 
can fail, although there is always existence of a unique solution which becomes 
instantaneously L°°. It also turns out that (USD can be interpreted as the 
gradient flow (as in [01 lAGSj l in the space of probability measures equipped 
with the 2-Wasserstein distance, of the energy functional 

(1.9) $(/x) = y |V/i|2 h = -A-V, 

which is also the mean field limit of the usual Ginzburg-Landau energy. The 
equation ()1.8p was also studied with that point of view in [AS] in the bounded 
domain setting (where the possible entrance and exit of mass creates difficul¬ 
ties). This energy point of view allows one to envision a possible (and so far 
unsuccessful) energetic proof of the convergence of (|1.2I) based on the scheme 
of Gamma-convergence of gradient flows, as described in |Se2j . 

A proof of the convergence of (HI), ([Ll or (11.31) to these limiting equations 
(Euler or (ll.SI) ) has remained elusive for a long time. 

The time-independent analogue of this result, i.e. the convergence of solutions 
to the static Ginzburg-Landau equations to the time-independent version of 
dH) - a suitable weak formulation of /zV/i = 0 - in the regime ^ 1 (the 
notation means -|-oo as e —)• 0), and without any assumption on the sign 

of the vorticity, was obtained in |SS2j . As is standard, due to its translation- 
invariance, the stationary Ginzburg-Landau equation can be rewritten as a 
conservation law, namely that the so-called stress-energy tensor is divergence 
free. The method of the proof then consists in passing to the limit in that 
relation, taking advantage of a good control of the size of the set occupied by 
vortices. This approach seems to fail to extend to the dynamical setting for 
lack of extension of this good control. 

On the other hand, the proof of convergence of the dynamic equations in the 
case of bounded number of vortices is usually based on examining the expression 
for the time-derivative of the energy density and identifying limits for each term 
(for a quick description, one may refer to the introduction of |ST2] L This proof 
also seems very hard to extend to the situation of Ag 3> 1 for the following 
reasons : 

• It relies on estimates on the evolution of the “energy-excess” D, where 
terms controlled in y/D instead of D arise (and these are not amenable 
to the use of Gronwall’s lemma). 

• Understanding the evolution of the energy density or excess seems to 
require controlling the speed of each individual vortex, which is diffi¬ 
cult when their number gets large while only averaged information is 
available. 

• The proof works until the first time of possible collision between vor¬ 
tices, which can in principle occur in very short time once the number 
of vortices blow up, so it seems that one needs good control and under¬ 
standing of the vortices’ mutual distances. 

Recently Spirn-Kurzke |KS2j and Jerrard-Spirn [.ISl] were however able to 
make the first breakthrough by accomplishing this program for dEZl) and dni 
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respectively, in the case where grows very slowly: <C (log log |log 

in the parabolic Ginzburg-Landau case, and ^ (log |login the two- 
dimensional Gross-Pitaevskii case, assuming some specific well-preparedness 
conditions on the initial data. Relying on their previous work |KSll |JSpl| , they 
showed that the method of proof for finite number of vortices can be made more 
quantitative and pushed beyond bounded N^, controlling the vortex distances 
and proving that their positions remain close to those of the points solving 
the ODE system (II.5p . and then finally passing to the limit for that system 
by applying classical “point vortex methods”, in the manner of, say, Schochet 
|Sch| . There is however very little hope for extending such an approach to larger 
values of N^. 

By a different proof method based on the evolution of a “modulated energy”, 
we will establish 

• the mean-field limit evolution of dEB-dlSI) in the regime |loge| <C 
Ne < 1/e 

• the mean-field limit evolution of (|1.2p in the regime 1 < 0(|loge|) 

both at the level of convergence of the supercurrents, and not just the vorticity. 
We note that the condition < 0(|loge|) allows to reach a physically very 
relevant regime: in the case of the equation with gauge, |loge| is the order of 
the number of vortices that are present just after one reached the so-called first 
critical field itself of the order of |loge| (cf. |SS5] i. 

Our method relies on the assumed regularity of the solution to the limiting 
equation, thus is restricted to limiting vorticities which are sufficiently regular. 
In contrast, although they concern very dilute limits, the results of [KS211.TSl] 
allow for general (possibly irregular) limiting vorticities. 

1.2. Our setting and results. 

1.2.1. Scaling of the equation. In order to obtain a nontrivial limiting evolution, 
the appropriate scaling of the equation to consider is 

(110) I + = + 

\ ri(-, 0) = u^. 

Here we have put both the 2D Gross-Pitaevskii and parabolic equations in the 
same framework. To obtain Gross-Pitaevskii, one should set a = 0 and /? = 1 
and to obtain the parabolic Ginzburg-Landau equation, one should set a = 1 
and /3 = 0. We will also comment later on the mixed case where one would have 
a and /3 nonzero, and say, + fP = 1, and we will describe the adaptation to 
the 3D Gross-Pitaevskii case as we go. In all the paper, when we write M”, we 
mean the whole Euclidean space with n being either 2 or 3. 

1.2.2. Limiting equation. Throughout the paper, for vector fields X in the 
plane, we use the notation 

X^ = (-X2,Xi) V^ = (-d2,di) 

and this way 

curlX = 81 X 2 — 82 X 1 = — divX"*-. 
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In the 2D Gross-Pitaevskii case, or in the parabolic case with <C |loge|, then 
the limiting equation will be the incompressible equation 


( 1 . 11 ) 


1 5fV = 2/3v-*- curl v — 2 q!V curl v -|- Vp 

in 

\ divv = 0 

in M^, 


with p the pressure associated to the divergence-free condition. In the parabolic 
case with comparable to |loge|, then without loss of generality, we may 
assume that 


( 1 . 12 ) 


A := lim 
£—^0 


I logs 

Ne 


exists and is positive and hnite. In that case, the limiting equation will be 


(1.13) 


dtv = A V div V — 2v curl v in 


This is a particular case of the mixed flow equation 

(1.14) SjV = — Vdiv V + 2/3v*‘curlv — 2Q:vcurlv in 

a 

The incompressible Euler equation is typically written as 
( 9fV-|-V • Vv + Vp = 0 in M” 

] div V = 0 in M”, 


(1.15) 


where p is the pressure. But when divv = 0 and n = 2, one has the identity 
(1.16) V • Vv = div (v (g) v) = v*" curlv -|- -V|v|^, 


so when a = 0, (|l.lll) is exactly the 2D incompressible Euler equation (up to 
a time rescaling by a factor of 2), with the new pressure equal to the old one 
plus |vp. Existence, uniqueness, and regularity of its solutions are well-known 
since Volibner and Yudovich (one can refer to textbooks such as |MBl IChl| and 
references therein). 

In the three-dimensional Gross-Pitaevskii case, our limiting equation will be 
the time-rescaled incompressible Euler equation rewritten again as 


(1.17) 


j dtv = 2dw (v(g)v-i|vp/)-k Vp 

in 

\ divv = 0 

in M^. 


Taking the curl of (11.111) . one obtains 

dtuj = 2div ((av*" + /3v)a;) uj = curlv, divv = 0, 


so in the case a = 0, /3 = 1, we of course retrieve the vorticity form of the 
Euler equation, and in the case a = 1, /3 = 0 and <C jlogej, we see that we 
retrieve the Chapman-Rubinstein-Schatzman-E equation (|1.8p . But in the case 
0 < A < -|-oo the curl of (|1.13l) is not (11.81) . In fact the divergence of v is not 
zero, even if we assume it vanishes at initial time, and this affects the dynamics 
of the vortices. This phenomenon can be seen as a drift on the vortices created 
by the phase of ti, a “phase-vortex interaction”, as was also observed in [BOSlj . 
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The existence and uniqueness of global regular solutions to (11.1111 can be 
obtained exactly as in the case of the Euler equation, see for example [Chi]. It 
suffices to write the equation in vorticity form as dtuj = A{lo) ■ Va; where A[uj) 
is a Fourier operator of order —1 (we do not give details here). In contrast, the 
equation (ll.lhl) is new in the literature, and for the sake of completeness we 
present in an appendix (cf. Theorem ED a result of local existence and uniqueness 
of solutions to the general equation (|1.14p . which suffices for our purposes. 
More results, including global existence of such regular solutions are established 
in |Du] . 

We note that the transition from one limiting equation to the other happens 
in the parabolic case, in the regime where is proportional to |loge|. This 
is a natural regime, which corresponds to the situation where the number of 
points is of the same order as the “self-interaction” energy of each vortex. 


1.2.3. Method of proof. Our method of proof bypasses issues such as taking 
limits in the energy evolution or vorticity evolution relations and controlling 
vortex distances. Instead, it takes advantage of the regularity and stability of 
the solution to the limiting equation. More precisely, we introduce what can be 
called a “modulated Ginzburg-Landau energy” 

(1.18) £e{u,t) = ^f |Vu - iujVe:V(t)P + {1 - \u\^)'tp{t), 

I l£ 

where v(t) is a regular solution to the desired limiting equation, and 


(1.19) 




p(t) - |v(t)|2 

-|v(t)P 


in the Gross-Pitaevskii case 
in the parabolic case. 


This quantity is modelled on the Ginzburg-Landau energy, and controls the 1? 
distance between the supercurrent = {iUe,Vuf} normalized by iV^ and the 
limiting velocity held v, because in the regimes we consider, the term A^|(l — 
|up)i/^(t) is a small perturbation (that term will however play a role in algebraic 
cancellations). 

One of the difficulties in the proof is that the convergence of js/N^ to v is not 
strong in L^, in general, but rather weak in due to a concentration of an 
amount 7r|loge| of energy at each of the vortex points (this energy concentration 
can be seen as a defect measure in the convergence of js/N^ to v). In order to 
take that concentration into account, we need to subtract off of £s the constant 
quantity 7rA^e|loge|. In the regime where 3> |loge|, then 7rA^£|loge| = o{N^) 
and this quantity (or the concentration) happens to become negligible, which 
is what will make the proof in the Gross-Pitaevskii case much simpler and 
applicable to the three-dimensional setting as well, but of course restricted to 
that regime. 

The main point of the proof consists in differentiating in time, and show¬ 
ing that ^£i;{us{t),t) < C{£e{ue{t),t) — vrA^ellogel), which allows us to apply 
Gronwall’s lemma, and conclude that if £^ — vrA^gllogej is initially small, it re¬ 
mains so. The difficulty for this is to show that a control in C{£^ — vrAigllogel) 
instead of C^/£l, is possible, even though the terms involved initially appear to 
be of order This is made possible by a series of algebraic simplifications 
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in the computations that reveal only quadratic terms. An important insight 
leading to these simplifications is that one should work, by analogy with the 
gauged Ginzburg-Landau model, with gauged quantities, where N^v plays the 
role of a space gauge vector-field, and N^4> plays the role of a temporal gauge, 
(f> being defined by 


( 1 . 20 ) 




p(f) in cases leading to (II.lip . ()1.17p 

Adivv(f) in cases leading to (11.1311 . 


The idea of proving convergence via a Gronwall argument on the modu¬ 
lated energy, while assuming and using the regularity of the limiting solution 
is similar to the relative entropy method for establishing (the stability of) hy¬ 
drodynamic limits, first introduced in |Yau] and used for quantum many body 
problems, mean-held theory and semiclassical limits, one example of the latter 
arising precisely for the limit of the Gross-Pitaevskii equation in |LZ2j : or Bre- 
nier’s modulated entropy method to establish kinetic to fluid limits such as the 
derivation of the Euler equation from the Boltzmann or Vlasov equations (see 
for instance |SRj and references therein). 

In the point of view of that method, |Vu — mV^vp is the modulated 
energy, while \dtu — iuN^(p\‘^ is the modulated energy-dissipation. 


1.2.4. Well-posedness of the Cauchy problem. The equations dEU) and dEl are 
shown in [G^IG^ to be globally well-posed in the natural energy space 

(1.21) E = {ue e L‘^(R^),\u\^ - 1 £ ^^(M’^)}. 


This is the setting we will consider in dimension 3 and corresponds to solutions 
which have zero total degree at inhnity. But in general this is too restrictive 
for our purposes: the problem is that, when working in the whole space, the 
natural energy is infinite as soon as the total degree of u at infinity is not zero. 
It thus needs to be renormalized by substracting off the energy of some fixed 
map which behaves at infinity like u^, i.e. typically like for example. 
To be more specific, in the two-dimensional case we consider as in [BSl IMi| . for 
each integer D, a reference map Ud, which is smooth in and such that 

/ \ D 


( 1 . 22 ) 


outside of B(0,1). 


The well-posedness of the Cauchy problem in that context was established 
in [BSj in the Gross-Pitaevskii case and |Mi] in the mixed flow (hence possibly 
parabolic) case : they show that given (in fact they even 

consider a slightly wider class than this), there exists a unique global solution to 
(jl.lOp such that Us{t) — Ud^ £ and satisfying other properties 

that will be listed in Section EH This is the set-up that we will use, as is also 
done in |JSp2| . Without loss of generality, we may assume that > 0. In 
the Gross-Pitaevskii case, we will allow for (the total degree) to be possibly 
different from (the total number of vortices), which corresponds to a vorticity 
which does not have a distinguished sign. For simplicity, we will then assume 
that D^/N^ = d < 1, a number independent of e. In the parabolic case, we 
need to have a distinguished sign and we will take = N^. 
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Let US point out that the use of the modulated energy will simplify the proofs 
in that respect, in the sense that it naturally provides a finite energy and thus 
replaces having to “renormalize” the infinite energy as in [B.TSl IMil |JSp2| . 

1.2.5. Main results. We may now state our main results, starting with the 
Gross-Pitaevskii case. In all the paper, we denote by the space of func¬ 
tions which are bounded and Holder continuous with exponent 7 , and by 
the space of functions which are bounded and whose derivative is bounded and 
. We use the standard notation a <C 6 to mean lima /6 = 0, and the standard 
o notation, all asymptotics being taken as e — 0 . 

Theorem 1 (Gross-Pitaevskii case). Assume satisfies 

(1.23) |loge| <C <C - as e ^ 0. 

Let {a£}e>o 6 e solutions to 

f iN^dtUs = Auc + ^(1 - m M" 

(1-24) < 

\ •«£(•, 0) = u°. 

If n = 2 we assume € Ud^ + LI^(IK^) where we take 0 < 
with D^/N^ = d, Uds is as in (ll.22j> . and y is a solution to (ll.lljl . such 
that v{0) — d{VUifiUi) € L^(M^), v(t) S L°°(M_|_, curlv(t) G 

L~(M+,L 1 (]R 2 )), 

If n = 3 we assume £ E as in ()1.21h with Au^ G L^(M^), and V is 

a solution to (jl.l7p such that v(t) G L°°([ 0 , T], C°4(]^3^ p g 

L°°([0,r],L°°(M3)nL2(M3)) and curlv(t) G L°°([ 0 ,T],^^(M^)). 

Letting he defined from v{t) via (11.181) . assume that 

(1.25) f,(aO,0) <o(iV2). 

Then, for every t > 0 (resp. t < T), we have £fiue{t),t) < o{N^), and in 
particular we have 

(1.26) — :=- — -^ V strongly in ). 

If we know in addition that is bounded in L°°(EE)) uniformly in e, 

then the convergence is strong in L^(M”). 

The restriction S> |loge| is a technical obstruction caused by the difficulty 
in controlling the velocity of the individual vortices because of the lack of control 
of /jjjn IdtUsl"^. On the other hand, the restriction <C 4 seems quite natural, 
since when is larger, the modulus of u should enter the limiting equation, 
giving rise to compressible Euler equations. On that aspect we refer to the 
survey [BDGSSj and results quoted therein. 

In the parabolic case, we have the following result 

Theorem 2 (Parabolic case). Assume satisfies 

(1.27) 1 <C A^e < 0(|loge|) as e —0, 
and let A be as in dLH. 

Let {a£}e>o be solutions to (ll.lOp with 13 = 0 and a = 1, associated to initial 
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conditions G where is as in (|1.22p . Assume v is a solution 

to (11.111) if <C |loge|, and to (I1.13p otherwise, such that v{0) — {VUi,iUi) G 
curlv(O) > 0, and belonging to L°°([0, T], for some 7 > 0 

and some T > 0 (possibly infinite). Letting Eg, be defined from v(i) via p.l 8 l) . 
assume that 

(1.28) < vriV,|loge| +o(A^2^. 

Then, for every t G [0,T], we have S,;{ue{t),t) < vrA^gllogel + o{N^) and 

(1.29) ^ ^ strongly in for p <2, 

and weakly in L^(M^) if in addition A < 00 , 

("1 o^^ (ViXgjZIXg) ; ■ rl /'Tn)2\ 

(1.30) — :=- — -^ V strongly in ). 

If we know in addition that is bounded in L“^(M+, L°°(M^)) uniformly in e, 
then the convergence of je/N^ is in the same sense as in (ll.29p . 

We note that in both theorems, we obtain the convergence of the solutions 
of (ll.lOp at the level of their supercurrents js , which is obviously stronger than 
the convergence of the vorticity He/Ate = curljg/iVg to curlv, which it implies. 

The additional condition on the uniform boundedness of Ue is easy to verify 
in the parabolic case: for example if the initial data satisfies |Ug| < 1, then this 
is preserved along the flow of (| 1 . 2 p by the maximum principle. 

The conditions (|1.25l) and (ll.28p are well-preparedness conditions. It is fairly 
standard to check that one can build configurations that satisfy them, for 
example proceeding as in |SS5[ Section 7.3]. 

In the parabolic case, for T 77^ (M^) and (I1.28P to hold, the initial 

configuration should have most of its vortices of degree 1, and thus curlv(O) 
must be nonnegative (it is automatically of mass 27r by the condition v( 0 ) — 
{VUi,iUi) G L2(r 2) so the assumption curlv(O) > 0 is redundant). We take 
advantage of these well-preparedness conditions as well as of the regularity of the 
solutions to the limiting equations in crucial ways. Since regularity propagates 
in time in these limiting equations, then the regularity assumptions we have 
placed really amount to just another assumption on the initial data. It is of 
course significantly more challenging and an open problem to prove convergence 
without such assumptions, in particular without knowing in the parabolic case 
that the initial limiting vorticity has a sign. 

The reason for the restriction Ne, < 0(|loge|) will become clear in the course 
of the proof: the factor of growth of the modulated energy in Gronwall’s lemma 
is bounded by CN^/\[oge\ and thus becomes too large otherwise. We are not 
even sure whether the formal analogue of (|1.13l) . i.e. the equation with A = 0 
(shown to be locally well-posed in [Du] i. is the correct limiting equation. 


1.3. Other settings. 


























MEAN FIELD LIMITS OF GINZBURG-LANDAU 


11 


1.3.1. The mixed flow case. With our method of proof, we can prove that if 

a > 0, /3 > 0, and = 1, the same results as Theorem [2] hold, with 

convergence to the limiting equation (ll.ll|l . respectively (I1.14p . under the ad¬ 
ditional condition 3> log |loge|, cf. Remark I4.1UI For a sketch of the proof 
and quantities to use, one can refer to Appendix C, setting the gauge fields to 0. 

1.3.2. The torus. We have chosen to work in the whole plane or space, but one 
can easily check that the proof works with no change in the case of a torus. The 
proof is of course easier since there is no need for controlling infinity, and there 
are no boundary terms in integrations by parts. However, this gives rise to 
a nontrivial situation only in the Gross-Pitaevskii case, since in the parabolic 
case we have to assume that the vorticity has a distinguished sign, while at 
the same time its integral over the torus vanishes. In the parabolic case, to 
have a nontrivial situation one should instead consider the case with gauge as 
described just below in Section 11.3.41 where the total vorticity over a period 
does not have to vanish. 

1.3.3. Bounded domains. On the contrary, working in a bounded domain en¬ 
tails significant difficulty in the parabolic case : one basically needs to control 
the change in energy due to the entrance and exit of vortices (one can see 
the occurence of this difficulty at the limiting equation level in |ASj 1. and the 
necessary tools are not yet available. 

1.3.4. The case with gauge. Our proof adapts well to the cause with gauge, 
which is the true physical model for superconductors, again in the setting of 
the infinite plane or the torus. The corresponding evolution equations are then 
the so-called Gorkov-Eliashberg equations 

I {dtu-iu<^>) = {Va)^u+ in m2 

\ — V$) = V-*-curl A-|-(in, V^tt) in M^. 

These were first derived by Schmid [Schj and Gorkov-Eliashberg [GEj . and the 
mixed flow equation was also suggested as a good model for the classical Hall 
effect by Dorsey [Do] and Kopnin et al. [KIK] . In this system the unknown 
functions are u, the complex-valued order parameter, A the gauge of the mag¬ 
netic field (a real-valued vector field), and $ the gauge of the electric field (a 
real-valued function). The notation denotes the covariant derivative V —zA. 
The magnetic field in the sample is deduced from A hy h = curl A, and the 
electric field hy E = —S^A-I-Vd*. Einally, cr > 0 is a real parameter, the charac¬ 
teristic relaxation time of the magnetic field, which may be taken to be 0. The 
dynamics of a finite number of vortices under such flows was established rigor¬ 
ously in jSpl] |Sp2]ISS41 IKSll[Til IST2| , in a manner analogous to that described 
in the case without gauge. A dynamical law for the limit of the vorticity was 
formally proposed in [GRSl lE^ . the analogue of (11.81) mentioned above. 

Natural physical quantities associated to this model are the gauge-invariant 
supercurrent 


je = AeUe) 
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the gauge-invariant vorticity 

/i£ = curl (je A^) 


and the electric field 

E, = -dtAe + Vd>,. 

In AppendixOwe explain how to adapt the computations made in the planar 
case withont gauge to the case with gauge, in order to derive the following 
limiting eqnations: if a = 0 or <C llogej, 


(1.32) 


(9tv — E = (—2av-|-2/3v-*-)(cnrl V-|-h)-|-Vp 
div V = 0 
—(tE = V -|- V-*“h 
, (9th = — curl E, 


and if a 7 ^ 0 and lime^o = A is positive and hnite, 

(9tv — E = {—2av + 2/3v-*-)(curl V -|- h) -|- Ay divv 


(1.33) 


-cjE = V -h V-^h 
(9th = — curlE. 


The corresponding results to Theorems [T] and [5] are then the convergences 


(1.34) ^ ^ cnrlv-Fh, 




N, 


N, 


curl A^ 
N, 




E 


in the case a = 0 and |loge| <C iV^ <C ^ to (ll.32p . and in the case (5 = 0 and 
1 <C A^e < 0(|loge|) to (|1.32p or (|1.33n according to the situation. 

We note that if (T / 0, (11.321) . resp. (11.331) . can be rewritten as a system of 
equations on only two unknowns v and h : 


(1.35) 


' 5tv -I- ^(v -|- V-*-h) = (— 2q:v -|- 2/3v-*-)( curl v + h) -|- Vp 
< div V = 0 
, adth = curlv-|-Ah, 


respectively 

f dtv -I- ^(v -|- V-*-h) = (—2av -|- 2/3v-*-)(curlv -|- h) -|- Ay divv 
(1.36) < 

[ adth = curlv-|-Ah. 

We can also note that taking the curl of (11.321) or (ll.33p gives (with fi the 
limiting vorticity as in (|1.34p ) 


(9tp = div ((2av-*- -|- 2/3v)/r) 

< adth — Ah + h = fj, 

divv = 0 or (9t divv = — A divv + div (—2 q!V + 2 /lv-’-)/i) — ^Adivv. 

This is a transport equation for p, conpled with a linear heat equation for h. 
As before, when choosing (5 = 0, this equation coincides with the model of 
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[CRSl IE2j when <C |loge|, but not if is larger and A is finite. 

The rest of the paper is organized as follows: We start with some prelimi¬ 
naries in which we recall some properties of the solutions and a priori bounds, 
introduce the basic quantities like the stress-energy tensor, the velocity and the 
modulated energy, and present explicit computations on them. 

Then we present the proofs in increasing order of complexity: we start with 
the proof of Theorem [T] which is remarkably short. We then present the proof of 
Theorem [21 which requires using all the by now standard tools of vortex anal¬ 
ysis techniques : vortex-balls constructions, Jacobian estimates, and product 
estimate. An appendix is devoted to the proof of the short-time existence and 
uniqueness of solutions to (ll.ljp . and another one to the computations in the 
gauge case. 

Acknowledgments: I would like to warmly thank Jean-Yves Chemin for his 
guidance with the proof of Theorem [3l and to thank Didier Smets for pointing 
out that the proof for the Gross-Pitaevskii case should work in dimension 3. I 
also thank Matthias Kurzke, Mitia Duerinckx, and the anonymous referees for 
their careful reading and very useful comments. 

2. Preliminaries 

In these preliminaries, we will work in the general setting of (jl.lOp (or its 
three-dimensional version) with a, /3, nonnegative satisfying = 1, which 

allows us to treat the Gross-Pitaevskii and parabolic (and mixed) cases at once. 
We note that with the choice ()1.20p . the limiting equations (jl.lip . ()1.13p can 
always be rewritten 

(2.1) 5tv = -|-2/Iv*-curl V — 2av curl V. 

Here and in all the paper, C will denote some positive constant independent of 
e, but which may depend on the various bounds on v. 

Also C'“^'’''^(M"') denotes functions that are a sum of derivatives of C''’'(M”) 
functions and C°°(M”) bounded functions, and C^''''’'(M”) is the same as C^’'>'(ffi"'), 
i.e. functions which are bounded and whose derivative is bounded and C'^. 

2.1. A priori bounds. 

2.1.1. A priori estimates on v. We first gather a few facts about the solutions 
V to (II.lip . (I1.13P and (I1.17P that we consider. 

Lemma 2.1. Let v be a solution to (jl.lll) in L°°([0, oo], C^’^(M^)) such that 
v(0) — d{'VUi,iUi) G L^(M^) and curlv G L°°(]R_(., L^(M^)). There exists a p 
such that (jl.lip holds and such that for any 0 < T < oo, v — diyUi,iUi) G 
L°°([0,T],L2(M2)), V G L°°([0,T],L4(m 2)), dtv G L°°([0, T], L2(m 2) nL°°(M2))^ 
p G L°°([0,T],C°’1(M2) nL2(M2)), and dfP,Vp G L°°([0, T], ^^(M^)). 

Let Y he a solution to (fTlTP such that v G L°^{[0,T],C^d^) n L^{R^)), 
dtv G L°°([0,r],L°°(M3) nL2(M3)) and curlv(t) G L°^{[0,T], L^iR^)). There 
exists a p such that (jl.l7l) holds and such that for any 0 < T < oo, p G 
L°°([0,T],C°’HM3)nL2(M3)) and dtp, Vp G L°°([0, T], ^^(M^)). 
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Let V be a solution to (|1.13p in L°°([0, T], 0 < T < oo, such that 

v(0) - {VUi,iUi) G L2(m2), We have v - {VUi,iUi) G L°°([ 0 , T], L2(m2))^ 
V G L°°([0,T],L4(M2)), divv G L°°([0,T],L2(M2)nL°°(M2)) and dtv,V divv € 
L2([0, T], L2(M2) n L°°(]R2)). Also for any t ^ [OjT], ^curlv(i) is a probability 
measure. 

Proof. Let us start with the case of p.ll|) . We first observe that such a solution 
exists, and that it corresponds to the solution belonging to the space -E 27 rd 
in the notation of [Chll Definition 1.3.3]. It is also known (cf. m) that 
the solution remains such that v(t) — d-^ G L°°([0, T], L2(m2^^(|o, 1))) and 

thus v(t) — (VC/i,iI7i) G L°°([0 , T], L2(m 2) (the case of general a and /3 works 
similarly). 

Since {VUi,iUi) decays like l/|x|, by boundedness of v and the Lp' character 
of V — d{'S/Ui,iUi), we deduce that v G L°°(M+, L‘^(M2)). 

The integrability of p is deduced from that of |v|2 by using the formula 

(2.2) Ap = 2/3 ^ div div (v (g) v) — + 2a div (div (v (g) v))-*- 

which means that Ap is a second order derivative of v (g) v. Since v G L^(M2) n 
( 70 , 1 ( 1 ^ 2 )^ this allows us to pick a pressure p such that p G L°° n L'^ (cf. |Chl[ 
p.l3]) and it is also i^y the assumed regularity of v. We also note that 
vcurlv G L2(M^) n L°°(M2) uniformly in time by boundedness of v, curlv, and 
the fact that curlv is integrable uniformly in time. On the other hand, we may 
write 

dw = curlv = V^-A”^ div (—2/3vcurl v + 2av*- curlv) 

where A~^ is the convolution with —^log, and with the above remark, we 
may deduce that dw remains in L2 (m2 ) p uniformly in time, and the 

result for Vp follows by using the equation. The same result follows for dtp by 
applying (9tA“^ to ()2.2h . 

The same arguments apply to prove the results stated for ()1.17h . 

Let us now turn to (|1.13p . It is proved in Theorem [3] that v(t) — v(0) G 
L°°([0,T],L2(m2)), dw G L2([0,T],L2(m 2))^ and divv G L°°([0,T],L2(m 2)). It 
also follows immediately that v(t) — {VUi,iUi) G L°°([ 0 ,T],L 2 (]^ 2 ^^ 
uniform character of v follows just as in the Euler case above. The fact that 
Vdivv G L2([0, T], L2(m 2)) follows in view of ()1.13D . The fact that ^ curlv 
remains a probability measure is standard by integrating the equation. 

Next, differentiating (II.I3p we find that 

(2.3) 5t(Vdivv) = AA(Vdivv) — 2V div (vcurl v). 

It can be found in |Ch2[ Proposition 2.1] that if u solves on M2 equation 
dtu — Au = / on [0,T] with / G L°°([ 0 , T], C'“ 2 >T'(]^ 2 ^^ initial data uq G 
C°’^(]R2) then 

(2-4) lkllL°°([0,T],Cr(E2)) < Ct (||'Wo||cr(K2) + II/IIl°°([0,T],C-2+7(r2))) 

(it suffices to apply the result there with p = p = oo and s = —2 + 7 and notice 
that the Besov space is the same as or (7°’^). Applying this to (12.3p . 
since the right-hand side is L°°([0,T],C'“2+7(]^2^^ we obtain that V div v G 
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L°°([0,T],hence L°°([0,T],Inserting into (11.131) and using 
that curlv G Ci L°° yields that dfV G L°°([0, T], L°°(M^)). 

□ 


2.1.2. Estimates on the solutions to (jl.lOp . 

Lemma 2.2. Assume and v{t) satisfy the assumptions of Theorem{l\ or[^ 
Let = 1 in the three-dimensional case of Theorem [I] and in the 

case of Theorem\^ Then for any T > 0 we have dtUe G L^([0, T], and 

for any t G [0,T], 

V{Ue{t) -Ud,),1- Vu£(t) - iUe{t)NeY{t) 

all belong to £^{ue{t),t) is finite, and 

j,(t)-iV,v(t)G(Ti + L2)(M"). 

Proof. First let us justify that clttte(t) G T^([0, T], T^(M”)). In the two-dimensional 
Gross-Pitaevskii case, since we assume £ lloe A- then studying the 

equation for '■= Ue — Un^, we find in [BSl Prop. 3, Lemma 3] that remains 
in (by semi-group theory) and thus dtUs = dtWs G by 

the equation. 

For the two-dimensional parabolic (or mixed flow case), we assume G 
17n + and (by decay of the energy for := m — Un) it shown in IMil 

Theorem 1] that dtu^ G Li([ 0 , T], 

In the case of the three-dimensional solutions to (USD, it is shown in [Gel 
Sec. 3.3] that when Au^ G T^(M^) then the solution belongs for all time to 
X'^ := {u G L°^fE?),D‘^u G T^(M^)} and thus in view of the equation dtUe G 
L-(R,L2(M3)). 

Let us turn to the two-dimensional cases. Following exactly [BSl Proposition 
3] or [Mil Lemma A.6], we get that 

^ [ |V(Ue(f) - t/De)P + \Ueit) - + (1 - \Ueit)\'^f < Cis,t) 

2 Jir2 

where C{e,t) is finite and depends on e, t and Ud^- We thus obtain the 
character of the Hrst three items. 

By Lemma 12.11 and = dN^ (resp. Dg = N^), we have that N^Y{t) — 
{XUDeAUnf) £ T^(R^). We may then write 

\XUe - iUeN^Y\ < |V(Ue - Ud^)\ + Ne\Y\\u^ - UdJ + \'^Ud^ - illDed^e^l- 

The first two expressions in the right-hand side are in L^(R^) by what pre¬ 
cedes and the boundedness of v. For the third quantity, we have that \VUDg — 
iUD^X^vl = \{XUDs,iUDs) — dLsv\ outside ofB(0,1) and is bounded inB(0,1), 
by definition of Ud, hence is in by Lemma 12.11 We conclude by Lemma 12.11 
that this term is also in T^(R^). The finiteness of 8e{ue{t),t) is then an imme¬ 
diate consequence of what precedes, the fact that v G from Lemma [2.11 and 
the Cauchy-Schwarz inequality. 

Lastly, for — N^y we write 

\je - N,y\ < \{VUD,,iUD,) - N,y\ + \V{Ud, - u,)||n,| + |V17 dJ|u, - UdJ 
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and we conclude by the previous observations, writing \ue\ = 1 + (lugl — 1) and 
using that (1 — \ue\)'^ < (1 “ that — N^v £ + L^(M"’). 

□ 

2.1.3. Coerciveness of the modulated energy. We check that the modulated en¬ 
ergy does control the quantities we are interested in. We have 

Lemma 2.3. The functional £^ being as in (jl.lSp and ip as in (I1.19p . we have 
for any Ue, R > I, 1 < p < oo, 

(2.5) \Vus-iu^Nev\‘^+ <£^{ue) + e‘^N^\\'ip\\l2, 


(2.6) [ \je - iVeVl < Cr^pWVUs - iUeN^v\\Rp(^R^) 

JBr 

+ Cr^p{£^N‘^\\'iP\\r2 + Ni,£\\y\\lco)£^{u^)^ + C££e{u^) + ||v||2,°o ||V'||^2 , 

and 

(2.7) / \je - N,V\^ < C{\\ue\\U + e^N^Ml^){£e{ue) + e^N^UWl,), 

JR" 

where C is universal, Cr depends only on R and Cr^p on R and p. 


In view of assumptions (ll.23p and (ll.27p and Lemma f2.11 the second term on 
the right-hand sides of ()2.5p and ()2.7p will always be bounded by o{N'^). Also, 
if llugllioo < C, the upper bound in (j2.7p is by C£e{ue) + o(l) and that in ()2.6p 
is by CR^p£e{Ue)^ + C££e{Ue). 


r I uuj. 




(1^ + _ |„.P) . 


2e2 

u 
1 
2 

For we write that 


Thus, using that f) £ Lf' hy Lemma l2.11 we have 

. IVUe - meiVeVpt)^ -h < ^^(ne) -F i . 

A ./pri 41:6 L /ron, 


( 2 . 8 ) 

|j£-iV^v| < \je-\UsfN^v\+Ne\l-\Us\‘^\\v\ = \{iUs,VUe-iUsN^v)\+Ns\l-\Ue\‘^ 

< \Ue\\VUe - iUsN^vl + iVe|l - |n£p||v| 

= \Vus - iUeN,;V\ + (luel - l)|Vne - iu^Nsvl + Ne\l - |n£p||v|. 

For (12.61) we integrate this relation over Br and use Holder’s inequality to get 
that for any 1 < p < oo, 


'Br 

+ ( [ (1-k.l)' 

'Br 


\je - A^Wl < Cr^p / IVUe - iUeNeV\P 


'Br 


'Bp 


\VUe - iUeNeVl'^ ) +Cr\\v\\l^N^ (I (1 - lUep)^ 


'Bp 
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and using that (1 — |«|)^ < (1 — \u\^)‘^ we are led to 

/ be - ^eV| < Cr^pWVUs - iUeNeY\\LP{BR) 

JBr 

+ Ce{£e{ue) + N,) {£,{ue) + N^UWl^Y 
hence (ESI) follows. The proof of E2D is a straightforward consequence of 

(ESI). □ 


2.2. Identities. In this section, we present important standard and less stan¬ 
dard identities that will be used throughout the paper. In all that follows v is 
a vector field, which implicitly depends on time and solves one of our limiting 
equations. 

2.2.1. Current and velocity. We recall that for a family {ueje, the supercurrent 
and vorticity (or Jacobian) are defined as 

je = {VUe,iUe) He = CUrlj^. 

Following |SS3j , we also define the velocity 

(2.9) Ve := 2{idtUe,VUe) 
and we have the identity 

(2.10) dtje = V{iUe,dtUe) + Ve- 

Taking the curl of this relation yields that dtHe = curll4. (In dimension 2, this 
means that the vorticity is transported by 1/^“, hence the name velocity). We 
also dehne the modulated vorticity 

(2.11) Hs := curl {{Vue - iUeN^Y, iue) + N^y) , 
and the modulated velocity 

( 2 . 12 ) Ve := 2{i{dtUe-iUeNs(p),VUe-iUeNsV) = Ve - N^vdt\Ue\'^ + N^(f)V\Ue\‘^, 
with (f) as in (|1.20l) . 

We will use the fact that for Ue solution of (11.1011 (resp. (|1.24p l we have the 
relation 

(y. 

(2.13) divje = Neiir. -r + i/3)dtUe,iUe), 

\\oge\ 

(resp. with a = 0 and /3 = 1) which is obtained by taking the inner product of 
(|1.10p or (I1.24P with iUe- 

2.2.2. Stress-energy tensor. We next introduce the stress-energy tensor associ¬ 
ated to a function u: it is the n x n tensor defined by 

(2.14) {Se{u))ki := {dkU,diu) - ^ (^|Vub + ^(1 - |ub)') 

A direct computation shows that if u is sufficiently regular, we have 

n 

div Se{u) := ^di{Se{u))ki = (Vu,Au-b —u{l - \uf)) 

1=1 ^ 
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SO if Ue solves (11.101) or (|1.24p . we have 

(2.15) dlY Se{Ue) = Ns{{ + i/3)dtUe,VUe). 

\loge\ 

We next introduce the modulated stress-energy tensor : 

(2.16) {Se{u))ki := {dkU - iuN^Vk,diu - iuNeVi) + N^{1 - |up)vfcV; 

- ^ ^|Vu - iuNev\^ + (1 - + ^(1 - ^ki, 

where 5ki is 1 if /c = Z and 0 otherwise. One can observe that 

(2.17) l'5e(u)| < |Vu — iuN^v\‘^ + ^^(1 — If “ 

and thus with Lemma 12.31 and the Cauchy-Schwarz inequality, we may write 

(2.18) / \S,{u)\<2Seiu) + N^\\l-\umMl,+2e^N^ml,. 

JR^ 

For simplicity, we will also denote for S^{u^) and Se for Se{ue), as well as 

(2.19) S'v := V (g) V — ^|vp/. 

Lemma 2.4. Let solve (|1.10l) or (ll.24p and v and (p be as above. Then we 
have 


( 2 . 20 ) 


div Se{Ue) = , {dtUe - iUeN^cj), - lUeN^w) + ^N^Ve - l3N^v{dtUe, Ue 

|loge| 2 

o N^CX I 19 TIT • / ^6^ I T \ 

-N^t. - r\u^\^V(j) + NeJeU -j-0 - dlV v) 


|loge| 


|loge| 


div 5v - iVe ((v • V)je -h {je ' V)v - V(j£ • v)) 


which in dimension n = 2 can be rewritten 


( 2 . 21 ) 

div5e(tte) = {dtUe - iu^Ne4>,'^Ue - iUeNeV) + ^N^Ve - l3N^v{dtUe, Ue:) 

|loge| 2 

-I- V*- curlv — Nej^ curlv — NeV^He 

+ N^v{divv- ^^^^\ue\^(j)) + Nejei .f^''^. (j) - divv). 

|loge| |loge| 

Proof. First, a direct computation yields 

Se{Ue) = Se{Ue) + - Ne{Y (g) je + je ® V - {je ■ v)I). 

Since we have the following relations for general vector fields v and j: 

(2.22) div (v (g) j) = j divv-|-(v • V)j, 


we deduce that 

( 2 . 23 ) dlY Se{Ue) = dlY Se{Ue) + dlY Sv 

- Ne {je div V V div je (v • V)je {je ■ V)v - V (je • v)) . 
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On the other hand, writing dtu^ = dtUe — iUsN^cp + iUeN^cp and Vn^ = Vtte — 
iUsN^v + iUeNfrW yields 


{dtUe, VUe) = {dtUe - iUeNe(p, VUe - iUeNeV) + N^jeCf} 

+ NeV{dtUe,iUe) - N^\Ue\'^V(j) 
and combining with (|2.9I) . (|2.13l) and (|2.15j) . we find 

Npa 


div Ss:{Us) = 


|loge| 


- iUeNe4>, VUe - iUeNe^) + NJs4> “ 

+ N^V divje - l3N^v{dtUe, U^) + 

Inserting into (j2.23h yields ()2.2np . In the two-dimensional case, we notice that 
we have the identities 


div 5v = V div v -|- v*" curl v 


and 

(v • V)j -|- (j • V)v — V(j • v) = curlv -|- v*- curlj, 
so (j2.21l) follows. 


□ 


2.2.3. Time derivative of the energy. Given a Lipschitz compactly supported 
function x(®)) and a sufficiently regular function 'if{x.,t) let us define 

(2.24) 

^e{u,t) = ^ X (^iVu - iniVev(f)p -h ^ + iVg (1 - . 

For simplicity we will most often write £e{ue) for £^{u^{t),t). 

Lemma 2.5. Let solve (|1.1UI) or (ll.24p and v satisfy the results of Lemma 
Then we have 

(2.25) 

^ie{Ue)=- [ + Vx ' (VUg - iu^N, dtuf) 

dt jRn I log e I 

+ / X (^£ V • dtY - Neje ■ dtY + {dtU^, iUe) div V - NeVe ■ v) 
JR" 

+ 11 xN^dt ((1 - \Uef){^P - |vp)) . 

^ JR" 

Proof. Since the solution is smooth and x is compactly supported, expanding 
the square, we may first rewrite as 


(2.26) 


£e{Ue) 



X 



2e2 J 





X^eje ■ V. 




















20 


SYLVIA SERFATY 


We then differentiate in time and obtain 
d f 1 

-T.^eiUe) = - xiidtUe, AUg + —Ue{l - • {VUe, dtUe) 

at £ 

+ [ x{N^^-dty-NJe-dtv-N^v-dtje) + \ [ X^'^dt ((1 - |wP)('0 - |vp)) . 

JKP- ^ JR" 

Inserting (|1.10l) and (|2.10l) and also writing Vn^ = Vug — iu^N^ ^ this 

becomes 

dti.e{Ue)=- / X |i ""“i \dtUe?‘ + Vx • ((ViXg - iUeN^W,dtUe) + N^w{dtUe,iUe)) 

jR" I log e I 

+ [ X (A^e V • dtV - Neje ' dtV - N^V • V {dtUe, iUe) - N^Ve ■ v) 
jR" 

+ l [ xN^dt ((1 - \Ue\^)i'lp - |v|2)) . 

^ jR" 

With an integration by parts, we find that two terms simplify and we obtain 
the result. □ 


3. The Gross-Pitaevskii case : proof of Theorem [T] 

In this section, we consider the Gross-Pitaevskii cases, in which v solves 
(II.lip with a = 0 and /3 = 1 or (I1.17p . and for which divv = 0, cj) = p and 
Tp = p— |vp. Below, we apply the result of Lemma 12. 5 1 with these choices. First, 
we insert the equation solved by v and the result of Lemma 12.41 to obtain the 
crucial step in our proof, where all the algebra combines. We note that in view 
of (ll.lbp . (II.lip and (I1.17P can both be written as 

(3.1) 9fV = 2 div 5v-|-Vp 
with the notation (I2.19p . 

Lemma 3.1. Let solve (ll.24h and v solve (II.lip or (I1.17p according to the 
dimension. Then we have 

(3.2) dtie{ue)= [ xi^Seiue) :Vv - N^{1 - \us\‘^)dt{\v\‘^ 

jR" ^ 

- [ '^X- (C^Ue - iU6Ney,dtUe) + Ns{NsV - je)p -2Sev), 
jR" ^ J 


where for two 2x2 matrices A and B, A: B denotes '^ki^kiBki- 

Proof. Starting from the result of Lemma 12.51 the first step is to insert (13.ip , 
which yields 

(3.3) [ X (A^e V • dtv - Neje ' ^jv) = [ x^dNeV - je) ' (2 div S’v -k Vp). 

jR" jR" 
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In order to transform the (first order in the error) term {N^v — j^) ■ div into 
a quadratic term, we multiply the result of Lemma 12.41 bv 2xv, which yields 

(3.4) -f xiVeK • V =-2 / xdiv5£(u£) • v +2 / xiV|div5v-v 

JR” JR” JR” 

- 2iV£ / X ((v • V)j£ + (je • V)v - V(j£ • v)) • V - [ xNeH'^dt\us\‘^ . 
JR” JR” 

But by direct computation, one may check that for any vector fields v and j we 
have 


(3.5) ((v • V)j + {j • V)v - V(j • v)) • V = -divSv -i 

and applying to v and v, we easily deduce that 

div 5v • V = 0. 


Therefore, inserting (13.5p applied with v and j^, (13.3p and (13.41) into the result 
of Lemma 12.51 and noticing several cancellations, we obtain 

d f 

-T^^eiUe) =- Vx • (VUe - iu^N^V, dtUe) 
at 

+ [ xiNsiNeV - je) 2 div Si;{Us) -v) 

JR” 

+ (iVs'(|vp5i(l - |u,|2) + ^N^dt ((1 - |u,l2)(V. - lv|2))) . 

Integrating by parts and using (|2.13p . we have 

/ X^eiNeV - js) • Vp = - [ NeVx ' (A^eV - js)p + [ N^x{dtUs,Ue)p. 

JR” JR” JR” 

Inserting into the previous relation and collecting terms, we are led to 
d f 

-TL^eiUe) = - Vx • ((VUg - iUeNeV, dtUe) + Ne{N^V - js)p) 
at J^n 

-2 xdiv^£(u£)-v 

JR” 

+ XNe (^5t(l - \Ue\^) (^|v|2 - ip^ + iSf ((1 - \Us\^)((llt - |v|2))^ . 

Since we have chosen i/; = p — |vp we see that the terms involving 9^(1 — \us'\^) 
cancel and we get the conclusion. □ 

We may check that all the terms in factor of x or Vx in (13.21) are in L^([0, t] x 
M”), thanks to (I2.17p . Lemmas 12.11 and 12.21 Inserting for x in (13.21) a sequence 
{Xk}k of functions bounded in C^piR") such that Xfc —^ 1 and ||Vxfc||L°°(M") —^ 0, 
we may thus integrate O in time and, using the finiteness of £c{uc{t),t) given 
by Lemma 12.21 let /c —)• oo to obtain by dominated convergence. 


- S,{ul0) = f [ 2~Se{Ue) : Vv - ^^(1 
JO JR” 


-Km(|v|2-0. 
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Inserting (12.181) . the bounds given by Lemma r2.ll and using the Cauchy-Schwarz 
inequality to control — \ue\^) hy in view of (12.5p . we arrive at 

£eiueit),t) - fe(u°,0) 

<C £eiu,) + CN^ ^ ^ eV£,{u,) + Ce^Ni^ 

<C f £s{u,) + Ct£^N^, 

Jo 

where C depends only on the bounds on v. Applying Gronwall’s lemma and 
using (II. 28^ . we finally obtain that 

£,{u,{t).t) <Ct{£e{ulQ) + o{Nl)) 

which is o{N^) with the initial data assumption (ll.25p . The conclusions of the 
theorem follow in view of either (12.61) or (j2.7l) . 


4 . The parabolic case 


In this section, we turn to the proof of Theorem [21 which only concerns the 
dimension n = 2. Throughout we assume that a = 1, /3 = 0 and (I1.27P holds. 
In the rest of the paper, we let XR be 


( 4 . 1 ) 


Xr{x) = < 


log|3;|-logR2 r> ^ ^ d2 

logR-logR^ tor K<\X\<K 

Xr{x) = 1 for |x| < R 
Xr{x) = 0 for |x| > R^. 


With this choice we have 


( 4 . 2 ) ||Vy;_R||i;,oo(R 2 ) ^ 0 ||Vxr||l 2 (]r 2 ) —0 as R^oo. 


4.1. A priori bound on the velocity. We dehne as the maximum time 
t < min(l,T) (where T > 0 is the time of existence of the solution to the 
limiting equation) such that 

(4.3) £eiu^{t)) < vrAellogel + for all t < T^. 

Our goal is to show that £e{ue{t)) < 7rAe|loge| + o(A|) for all t < T^, which 
will imply that Tg = min(l, T). 

Let us start with a crude a priori bound on the time derivative of u^- 

Lemma 4.1. Assume (jl.27l) and £s{u^) < nAgllogel + o{N^). Assume v 
satisfies the results of Lemma \2.1[ Then 

r [ \dtu,\^<CN^\loge\^ 

Jo JR2 

where C depends only on the bounds on v. 
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Proof. Let us return to the result of Lemma 12.51 with the choice = — |vp and 
X = XR as in (BH). By (j2.9p ^ AV0 li3iV6 ^ \/^ \/ 

iUeN^y) -\-2N^Y{dtUe^ and inserting into the result of Lemma [231 we obtain 

dtSeiUe) = - [ XR n^^ I \dtUe\‘^ - f Vxr ' C^Ue - iUeNsV, BtUe) 

JlR2 I log e I JlR2 

[ XR (-iv|(l - \Ue\‘^)dt\v\‘^ + Ne{dtUe,iUe)divv 
JK 2 

+Ne{NsV - je) ■ dtV + 2Ne{idtUe,VUe - iUeNeV) ■ v) . 


+ 


We next observe that again all the terms in factor of XR or Vxi? are in L^(M^) 
thanks to Lemmas 12.11 and 12.21 We may then integrate in time and let i? —oo 
to obtain 
cTe 


ri 


N,, 


-\dtUe\^ 


R 2 |loge| 

= £e{u^^) - £e{ue{T^)) + [ [ N^idtUs, iu^) div V + Ns{N^Y - j^) ■ dtY 

Jo Jr'^ 

+ 2Ne{idtUe, VUs - iUsNsv) ■ V + o(l), 

where we also used (ll.27p and (14.3h to control all the terms containing (1 —lu^p). 
Next we insert (12.8D to obtain 


n 

Jo 




’0 JR^ |loge 

fT, 


■\dtUe\^ < £e{ul) 


+ f [ fACe|9tUel| divv| + |1 - lUelllStUell divv| + iVgIVUe - iUeiVeVllcliVi 

Jo 4k2 V 

+ iVell - \Ue\\\VUe - iUeN^v\\dtv\ + N^\l - |u£p 11 v| |9tv| 

+ 2N^\dtU^\\VUe - meiVev||v|^ + o(l). 

Using |1 — lugll < |1 — lugpl, the Cauchy-Schwarz inequality, (11.281) . (14.3p . the 
L°° n character of div v and dtv given by Lemma 12.11 the boundedness of v 
and Lemma 12.31 we deduce that 


/7 


N, 


-\dtUef 


K 2 |loge| 

< vriVellogel +C'iV£ [ {I + \\dtUe\\L 2 ){l + £eiue)) dt + o{N^). 

Jo 

Using (14.31) . bounding by 1 and using (ll.27p . we easily deduce the result. □ 


4.2. Preliminaries: ball construction and product estimate. The proof 
in the parabolic case is more involved than in the Gross-Pitaevskii case, and 
in particular it requires all the machinery to study vortices which has been 
developed over the years. Indeed, as explained in the introduction, we will 
need to subtract off the (now leading order) contribution of the vortices to 
the energy. This will be done via the ball-construction method (introduced 
in [HillJI]) coupled with the “Jacobian estimates” [JSlj (with precursors in 
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[BEll ISSlj ). Here we will need a lower bound with smaller errors, as in |SS5[ 
Theorems 4.1 and 6.1], coupled to an improvement due to |ST1| . 


4.2.1. Vorticity to modulated vorticity. Before doing so, we will need the fol¬ 
lowing result which connects the vorticity to the modulated vorticity. 

Lemma 4.2. Assume v satisfies the results of Lemma \2.1\ and is such that 
£e{ue) < oo. Let = curlje and fie = curl((Vii£ — iUsN^vfiUe) + N^v) as in 
(|2.11l) . We have that G with 

(4.4) l/lgl < 2|Vti£ — itieA^gVp + |1 — l^^ePlI curlvj 
and 

(4.5) f jle = 2 TTNe. 

JR2 

Moreover, for any ^ G iL^(M^), we have 


(4.6) 


/ Cihe-he) 

./r2 


< CWVaLHRfiNfiey/i;^ + Ce^N^). 


Proof. First, a direct computation gives that ( 1131 ) holds, and it follows imme¬ 
diately with Lemmas 12.11 and 12.21 that G L^(]R^). We may then write, with 
XR as in gH) 


/ fL^= lim / XRhe 

2 R 2 R-^OO 

= - lim / V-^XR-{{'^Ue-iUeNeV,iUe) + Ns{'v-{'^Ui,iUi) + {VUi,iUi)). 
Jr2 

Since — iu^N^v and v — {VUifiUi) are in Lfi by Lemmas 12.11 and 12.21 the 
corresponding terms tend to 0 as ii —>• oo. There remains 



N,V^XR-{^Ui,iUi) 

= lim [ NsXRCUi\{'VUi,iUi) = 2TrNe 
2r2 


by choice of Ui. For (|4.6I) we observe that by a direct computation, we have 
fie — he = curl ((1 — \ue\‘^)Nev). Thus, for any ^ G i/^(M^) we have 



he) 



'W£|^)V-^C • V 


< Ne\\VaLHR^){Ce^£e{ue) + e^Ni), 


where we used Lemma 12.31 


□ 


4.2.2. Jacobian estimate for unbounded domains. The next lemma is an infinite 
domain version of the estimate of [SS51 Theorem 6.1]. For this lemma, we 
temporarily use the notation /Xg with a slightly different meaning. 
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Lemma 4.3. Let be an open subset o/M^, and let = {x G Q, dist (x, dfl) > 
e}. Let tie : 0 —)• C and : 0 —)• Assume that {Bi}i is a finite collection 

of disjoint closed balls of centers Oj and radii ri covering {Hugl — 1 ] > 5 } n , 
and let di = deg(iXe, (95^) if Bi C and di = 0 otherwise. Then, setting 

He = curl ((VUe - iAeUe, iUe) + Ae), 
we have, for any f, G C'c’^(n), 


(4.7) 


[ fine- diSafi 
Jn , 


- J^2\VUe-iUeAe\‘^ + \l-\Ue\‘^\\curlAe\ + 

where C is universal. 


2\2 


2e2 


Proof. As in [SS51 Chap. 6 ], we set x ^ ^ I^+ to be defined by 

x{x) = 2x if X G [0, |] 

X{x) = l ifxG[J,|] 

X(x) = 1+ 2(x - 3/2) ifxG[|,2] 

xix) = X if X > 2 . 

We then let Ve{x) = This is clearly well-defined, with |xe| = 1 outside 

of UiBi and from jSSSl Chap. 6 ] or direct calculations, we know that 

(4.8) IKVWe; - iAeUefiUe) - {VVe - iAeVe,iVe)\\Li{n) 

— 3||t “ I W| ||i,2(q/| VUg — iAeUe\\i,2(^fi'j, 

and 

(4.9) |1 - |Xg|| < |1 - Irigll, \VVe-iAeVe\<2\VUe-iAeUe\. 

Letting fie = curl ((Vug — fAgXg, ixg) -|- Ag), we have 


(4.10) 


iihe- fie) 


^ ((^^e iUeAe,iUe) (VXg fAgXg, ZXg)) 


< 3||V^||ioo(f^)||l — |ttg|^||j;^2(Q)||VUg — fAgUg||^2(Q). 

Next, we note that fie vanishes wherever |ug| = 1, and thus as soon as ||rig| —1| < 
Thus, by property of the balls, we have Supp fie r\ C UiBi (recall the 


2’ 


definition of in the statement of the lemma). We also have that whenever 
Bi C it holds that J^fie = 2TTdi (see 


Lemma 6.3]). Writing ^ as 
(,{ai) + O(rj)||^||c-o,i in each Bi, we conclude, exactly as in the proof of |SS5[ 
Theorem 6.1]. The only point that is a bit different is we need an analogue of 
Lemma 6.4] to bound J2i \He\ which works on an unbounded domain. 


For that, we may check by direct computation that 
(4.11) He = 2{VUe - iAeUe,i{VUe - lAeUe)) + (1 “ ]rtg^) curl Ag, 


so 


/ \He\ < 2 / 2lVUg - iAeUe 

i/ Q «/ Q 


-I- jl — |ug]^ll curl Ag 
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and the same holds for with a factor 4 in front, in view of (14.91) . We thns 
obtain that 



277'^didaj 
i 


< + e)||^||co,i / 2|Vn£ - + |1 - luepll curlAel 

j Ju,BiU{n\n^) 

+ C'IICIIc'0>l I|1 “ l^£PllL2(E2)l|V^£ - *^£^£||l2(K2), 


where C is universal, and we easily deduce the result. 


□ 


4.2.3. Ball construction lower bound, sharp version. In the result below, we use 
the Lorentz space L^’°° as in Eig, which can be defined by 

(4.12) 11/11^2,oo(k 2 ) = sup [ I/I 

|E|<oo JE 

where \E\ denotes the Lebesgue measure of E. 


Proposition 4.4. Assume £e{ue) < 7rA^e|loge| + where Ng, satisfies (ll.27p . 
and V E with curlv E Then there exists Eq such that for any 

e < Eq, the following holds. There exists a finite collection of disjoint closed 
balls {Bi = B{ai,ri)}i such that, letting di = deg{ue,dBi), the following holds 

(1) Hin < 

(2) {x, ||ue| - 1| > 4} C Ui5(aj,ri). 

( 3 ) 

\ f |VUe - iUeNsV\‘^ + ^ > TT ^ |dj| |log gj - o{NI). 

( 4 ) 


||Vn, 

( 5 ) 

(4.13) 


- m£A^£v|||2,oc(R2) < C - 7rY^\di\\logE\ + ^ |di|^^ 


Letting fie be as in ()2.1ip . for any ^ E C| 




, we have 




<o(l)llellc0.r. 


+ 0(iV2). 


Proof. The result of |SS5[ Theorem 4.1] applied to Ue, = A^^v and a = 3/4, 
provides for e small enough, for any e^/‘^<r<la collection of disjoint closed 
balls B{r) covering {x, ||u£| — 1| > such that the sum of the radii of the 

balls in the collection is r, and such that denoting D := YlseBir) I^bI, we have 


1 

2 


'UsgB(r) 


|VUe 


iUsN^v\‘^ + r'^\Ne curlvp + 


2^2 


(i-i^£n 

2e2 


> ttD 



One notes that the fact that we are in an unbounded domain does not create 
any problem. Indeed, since < oo, this implies (cf. e.g. m Lemma 2.3] 
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or IHLl Lemma 3.5l)that there exists a radius Re such that luel > i outside of 
B{<d,Re) and 

[ \VUe - iUeNeVl'^ + ^ (1 - < 1 

JB{ld,ReY 

i.e. the remaining energy is smaller than the desired error. This way the con¬ 
struction can be applied in B{0, Re + 1) only, yielding always a finite collection 
of balls covering {luej < |}. 

We first apply this result with the choice r = r' = to obtain a collection 
of balls {Bj} with centers a'-, radii r'- and degrees d'-, covering {xjHugl —1| > 
and satisfying r'- < The estimates of |SS5j also yield that 

V |d'I < -r^ / \Vue - iUeNev\^ + curlv^ + 

J -^3 

Using that curlv is bounded, the fact that the upper bound on 

£e{ue) and ()1.27p . we deduce from this relation that |d' | < CNe- 

We next apply the above result with the choice r = . This gives a 

collection of balls {Bi} = {B{ai,ri)} of degrees di, satisfying items 1 and 2 of 
the proposition and 


^ f |Vu£ - iUeNevl'^ + e '^^\Ne curlv|^ -b ^ 

2 Ju^Bi 2e2 

^ / e-UW \ 

It is part of the statements of [SS41 Theorem 4.1] that the family I3(r) is in¬ 
creasing in r, i.e. here that the Sj’s cover the B^s. By additivity of the degree, 

we thus have — Ylj\d'j\ — CNe, and using U < and the 

boundedness of curlv, the desired estimate of item 3 follows. 

For item 4, we use [sm Corollary 1.2] which yields that 


(4.14) llVUe - iUeNeV\\\2,o 



This is essentially a strengthened version of the result of item 3, in which the 
difference between the two sides of the inequality is shown to bound from above 
||VU£ - iUeNeY\\\2,^(^^2y 

Let us now turn to item 5, which is an adaptation to an infinite setting of 
the Jacobian estimates, for instance as in [SS51 Theorem 6.1]. The reason we 
needed a two-step construction above is that the total radius of the second set 
of balls, ^ (which have to be chosen this large so that they contain enough 

energy) is not very small compared to |loge| when Ne is not very large, and 
thus the Jacobian estimate applied directly on these large balls would give too 
large of an error. 
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Letting ^ be a smooth test-function, we may write 


= ^ d'C(a' 

* j i 

Since B' cBi ^'j ~ Sj Mjl ^ CN^, we may write 


(4.15) 




< 


iieiicwE4( E i‘'L 


j.BrcBi 


< 


ii«iic«(E''‘)E E i<*ji<ciieiic»..e“’'^v = c.(i)ii«iioo... 

i i J,B'cBi 


Applying Lemma 14.31 with = N^v on and it yields that 


[ ^(271^4-5^/ -/ie) 

JK2 j 

< Cc^/^ll^llc-o.i [ [Vug - mgA'gVp + jVg|l - lugpll curlvl + ^ . 

Jm? 2 e"' 

Using that curlv G the upper bound on £lg(ttg), and combining with 

(j4.15ll . we obtain the desired result for 7 = 1. The result for 7 < 1 follows by 
interpolation as in [JSI], using that Y.j \d'j\ < CN^ and lAel < CAgllogel 
hence d'j6a'. - fie\\{cO)* < CN^. □ 

4.2.4. Ball construction lower bound, localized version. We will also need a less 
precise but localizable version of the ball construction. This can be borrowed 
directly from [JlllSairSSlj . and combined with the Jacobian estimate of Lemma 
Mi so we omit the proof. 

Lemma 4.5. Under the assumptions of Proposition EZl there exists Eq such 
that for all £ < Eq, there exists a finite collection of disjoint closed halls {Bi\i = 
{B{ai,ri)ii such that, letting di = deg(?Xg, 55^), the following holds 

( 1 ) EiU <e-vra. 

( 2 ) 

Vi, ^ / |Vug - mgAgvl^> 7r|di||loge|(l - 0 ( 1 )). 

J Bi 

( 3 ) For any 0 < 7 < 1 and any ^ G we have 

/ 

Jr 2 y 

We emphasize that these balls are not necessarily the same as those obtained 
by Proposition 14.41 


^ ^ didai he 


<o(l)Pllco.. 
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4.2.5. Consequences on the energy excess. We next show how the energy excess 
— "TiVgllogel controls various quantities, including the energy outside of the 
balls. 

Corollary 4.6. For any t < (where is as in ()4.3I) ). letting {Bi}i (de¬ 
pending on s,t) be the collection of balls given by Proposition 14-41 we have 


(4.16) 1 


\VUe - meiVeVp + 


( 1 - 


Uf 


| 2^2 


2^2 


< Seiue{t),t) -vriVgllogel +o(iV|), 


(4.17) 


for £ small enough 


N,<Y,\d^\<CN,, 

i 


(4.18) ||VUe - iU£A^£v||i2.oo(]R2) < CNe, 

and for any nonnegative ^ € C*^’'^(M^), 


< IICI|l°° {£e{Ue{t),t) - vriVellogel) + o(iV£) ||C||co.7 ■ 

Proof. First, applying item 5 of Proposition 14.41 with ^ = XR in (14.ip and 
letting R oo, we must have /]R2(2vr — fie) = Oe{C). Comparing with 

(USD, we deduce that Yhi di = Ne for e small enough. Subtracting the result 
of item 3 of Proposition 14.41 from £^{ue{t),t) and using Lemma 12.31 we then 
obtain (j4.16p . The upper bound in (I4.17P was proved in the course of the proof 
of Proposition 14.41 the lower bound is an obvious consequence of Yhi di = N^. 

The relation (j4.18p is a direct consequence of item 4 of Proposition 14.41 ()4.3p 
and (|4.17p . writing X). \di\‘^ < (J2i \di\f- 

Finally, for (I4.19P we use instead the balls given by Lemma [4. 5 1 From item 2 
of Lemma 14.51 we may write 


/R 2 \UiSi 


*? 5 /, 


\VUe - iUsN^v\‘^ + ^ 

2 , 


I Vue — iUsN^vr + 


2^2 


- 7r|di||loge|(l - o(l)) 


< 


£^{u^{t),t) -7r^|di||loge|(l + o(l)). 


Moreover, from the same argument as before we have di = and |di| < 
CN^ for these balls, hence since the terms of the above sum are all nonnegative. 
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adding to both sides 7rX^j(|di| — di)|loge|, we obtain 




\L 


2^2 


IVWe - iUeA^eVl^ + 


- iUeNM + 


2 e 2 


2 e 2 


- 7rdj|loge|(l - o(l)) 


< £eiue{t),t) - vriVellogeKl + o(l)). 


Next, separating the integral between UiBi and the complement, and using 
item 1 of Lemma 14.51 we deduce that 


(4-20) ]- [ ^|Vug - mgNgVp+ ^ < tt^ d^g(aj)|log£|(1 + o(l)) 

^ Jm2 ae ^ 

+ O(l)ll?llc0.r + ||?||l°° {£e{Ue{t),t) - 7rN£|loge|(l + o(l))) . 

Combining this with item 3 of Lemma 14.51 we deduce that (14.191) holds. □ 


4.2.6. Approximation ofv. We will need an approximation of v{t) based on the 
balls constructed via Proposition 14.41 


Lemma 4.7. Let v satisfy the assumptions of Theorem\^ and satisfy (|4.3p . 
For each t < Tg, letting {Bi}i be the collection of balls constructed in Proposi¬ 
tion 
that 


there exists a vector field v 


(depending on e and t) such 


(1) V is constant in each Bi, 

(2) for every G [Oj 1]; l|v —v||co,7 < — Ce~^^~"''>''^ where C 

depends only on v and 7 , 

(3) V — V has compact support. 


Proof. It is an adaptation of Proposition 9.6 of |SS5] . which we can apply to 
each component of v. We note that since the collection of balls is finite, we may 
replace v by v + x(v — v) where y is a smooth positive cut-off function which 
is equal to 1 on a large enough ball containing U^Bi and vanishes outside of a 
large enough ball. This makes v — v compactly supported without affecting the 
other properties. □ 


By continuity of (for fixed e) and of v, one may check that we may make 
v(t) measurable in t. While we have a good control on Vv, we have no control 
on dfV, and this is what prevents us from applying this method in the regime 
-Ng < 0(|loge|) in the Schrodinger case. 


4.2.7. Product estimate. Finally, to control the velocity of the vortices, we also 
need the following result, whose proof is postponed to Appendix 1X1 and which 
is an ^-quantitative version of the “product estimate” of |SS3j . 

We let Mg be a quantity such that 

(4.21) Vg > 0, lim = 0, lim = 0, lim = 0. 

e^o £^o M| £^o |loge| 

For example will do. In the statement below we do not aim 

at optimality, however we state the result in a way that would allow to go 
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beyond the range < 0(|loge|) that we are considering here (for example up 
to arbitrary powers of |loge| or quantities that satisfy the properties (|4.21l) '). 

Proposition 4.8. Let : [0, r] x —)• C. Let v be a solution to (jl.lip or 
(ll.bjp and (f) be as in (ll.‘2UD . Let be as in (|2.1‘2I) . Let X G C'‘’’^([0,t] 
be a spatial vector field. Set 

Pe = J (^j ^\dtUe - iUsN^fil'^+£fit)'^ dt, 

and assume < M^. Then for any A>1, we have, as e —>■ 0, 



where C depends only on the bounds on and v. 

The second line in the right-hand side is o(l), so is logM^/llogel. One can 
see that optimizing over A by taking 


_ / /q \dtUe iUsN^fil'^ \ 2 

\Io /r2 - iUeNeV) • X|2y 

yields a right-hand side in the form of a product (plus error terms), hence the 
name “product estimate”. 


4.3. Proof of Theorem [2l We now present the main proof. 


4.3.1. Evolution of the modulated energy. In the next result, we take as before 
= — |v|2 in the dehnition (12.241) . and insert the equation solved by v and 
(j2.21l) into (I2.25P to obtain the crucial computation. 

Lemma 4.9. Let Ug solve (11.101) and v solve (11.111) or (I1.13p . Assume that 
(jl.27p and (|4.3I) hold. Then, for any t < T,., we have 


(4.23) Efiufit), t) — £^(u^, 0) — Is Iv F Ie A Id F Id F Im F Iv F o(l) 
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where 


Is = 2 


Iv = - 


f! 

Jo JR' 

, 

Jo JR 2 

rt 


Se : Vv" 


NeVe • V 


Ie = - 


Id = - 


fi 

Jo JR 2 

/7 

Jo JR 2 


2Ne\'v\^fi,e 




\dtUe - iUeN^(j)Y - Ns{dtUe - iUsNsCj), iUe)idivv - 




0 Jr^ |loge| 

- 2 ,,^^ , {dtu^ - iUeNsCj), VUe - iUeNsv) ■ v-* 


|loge 


■ 0 ) 


|loge| 


Ih = 


Im — 


Iv = 




lO JR2 


2NsV-^ ■ {je - iV£v)(-- -4> - divv) 

|loge| 




2N^{je - NeV) • (v - v) curlv + 2N^{v - v) • vjle. 



0 JR 2 


Proof. Again, we start from the result of Lemma 12.51 applied with if = —|vp. 
The first step is to write 

(4.24) - f 

JR 2 I log £ I 

^ ~ f I {\dtUe - - N^\Ue\‘^(l)^ + 2N^(t){dtUe,iUe)) . 

Jr 2 I logs I 

The second step is to use (12.ip to write 


(4.25) f x[Ne^ ■ dt^ - Neje-dt^) = f x^£(“2|vpcurlv + v-V0) 

JR2 JR2 

- xNejs ■ {- 2 vcutIv+ V(j)). 

JR2 

On the other hand, integrating by parts and using (j2.1,3h . we have 


(4.26) / xNeiN,Y-j,)-V(P 

JR2 

= - [ AeVx • {NeV - je)(f> “ [ N^X ( divv - -7—iUg) ) (j). 

7r2 7r2 V |loge| J 


Next, we would like to transform the linear term JJg 2 x(~2A^|vp curl v + 
2N^j^ -vcurlv) into a quadratic term plus error terms. For that, we would like 
to multiply the result of Lemma 12.41 by 2v-’-, which after integration by parts 
leads to terms in : Vv*-. Using v rather than v leads instead to integrals 

that live only outside of the balls since Vv = 0 there by item 1 of Lemma 14.71 
These terms will thus be controlled by the energy outside of the balls, i.e. the 
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excess energy, as desired. This creates an additional set of error terms in v — v, 
which we will control thanks to item 2 in Lemma 113 

So as explained, let us multiply the result of Lemma 12.41 bv 2xv“’-, where v 
is given, for each time t, by the result of Lemma 14.71 This yields 


/ 


2 yv-^ • diy Ss{ue) 


L 


= X 2 




|loge 


■{dtU^ - iUeNe4>, VUe - iUeNeV) ] ■ V 


+ / X lv|^ curlv — 2A^£j£ • v) curlv — 2A^£|vp/r£) 

iK 2 

+ X {2N^ (v — v) • V curl v — 2• (v — v) curl v 
JR 2 


-2iVe(v - v) • V/ie) 


+ [ xfevv^(divv- 
iR2 V 


\Ue\‘^(j)) + 2NJs ■ V-^(t 


— divv) ) . 


|loge| |loge| 

Inserting (j2.12p , (14.241) , (j4.25p and (14.261) into the result of Lemma 12.51 applied 
with = — |vp, taking advantage of the cancellations and using one integration 
by parts, we obtain 

(4.27) 


—= - / X7j-f 

dt J ^2 I log e I 


- iUeNe<P\^ - N^\Ue\^\<P\^ 


- Vx • ((Vue - iueNeV, dtUe) + N^{N^y - js)(l> - 2Sey-^) 
4R2 

+ ^2 ^ (dXe{dtUe,iUe) div V - </> div V^ 

+ [ X (Ne{-Ve - Nevdtlu.l^ + N,cj)V\u,\^) • V + 25, : Vv^ 

JR2 ^ 


/ 

Jr 


+ X 2 


N, 


■{dtUe - iUeNe4>, Vue - iUeN^v) ■ - 2A^,|v|^//, 


|loge 

+ / X (v — v) • vcurlv — 2N^j^ • (v — v) curlv — 2A^,(v — v) • v/r,) 

JR 2 

+ f x 2 NeV-^ ■ ({js - iV,v)(—- divv) + —- |u,p)(/>v^ 

7r 2 V |loge| |loge| ) 

-[ xN!dt{{l-\us\^)\v\^). 

4R2 

Let us make three transformations to this expression. First, let us single out 
the terms 

[ xN^{-M^dt\u,\^ + • V|u,|2 - dt{{l - |u,|2)|vp) 

4R2 

= f x^e (1 - l'«£p)(-5t|v|^ + V((> • V + (()divv) + [ fV|Vx • v(l - |u,p)0 
7R2 Jr2 
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(where we have used an integration by parts). 

Second, let us replace {dtu^, iu^) by {dtUe — iUeN^cj), iue) + N^\us\‘^(l) in (I4.27h . 
this leads to a cancellation of the terms in (1 — |ttep)(/)div v and in 
Third, owing to (I4.6p let us replace —2 Jjg 2 X-^e|vp/Ue by 

-2 xNsM^fie + O (dl Vx||l2 + II Vv|U^ ||v||i2)iV,(ev/4(^ + e^N^)) 

and the same for —2 Jjg 2 — v) • v/ij. Both give rise to o(l) error terms by 
(fTTZll and M . 

After these substitutions, let us integrate in time and take x = XR- We may 
check via Lemmas 12.21 and EH and the fact that v — v is compactly supported, 
that all integrands in factor of XR or ^XR ^re in L^. Letting R oo, we then 
get the conclusion. □ 

We next turn to studying these terms one by one. We will show that Id,Im-, Iv 
are all negligible terms, while Iv,Ie and Id recombine algebraically thanks to 
the product estimate, to give a term bounded by the energy outside the balls, 
as does Is- 

4.3.2. The negligible terms. Let us start with Id- In the case <C |loge|, from 
(jl.20p we have 4> = p and we also have divv = 0. Then, in view of (|4.3H and 
(EH) we may bound 

\Id\ < 2 / Adpr^lvllje - A^vIIpI 
Jr2 I log e I 

_ /Y2 

< C'l|v||L°°||p||i, 2 VvrA^e|loge| = o{N^), 

where we used the boundedness of v hence of v and the character of p (see 
Lemma EH- la the case —)• A > 0 finite (see (11.121) 1. we have (j) = A divv 

and 

(4.28) 7 -— ^(j) — divv = (--—— — 1) divv = o(l)| div v|. 

|loge| |loge| 

We again conclude easily that |/d| < o{N^) in that regime too. 

The term 1^ is easily seen to be o{N^), using the Cauchy-Schwarz inequality, 
()4.3p and Lemma [2.31 and the integrability of dtv, p, Vp, divv and Vdivv 
provided by Lemma l2.ll 

To bound the first terms of R we first use ()2.8p and the fact that v € 
to get 

(4.29) 

/ {je - iVeV) • (v - v) curl V 
4r2 

+ |1 — iMglllVwg — meiVgvl + iVg |1 — |u£p|| curlv|^ 

Using the Cauchy-Schwarz inequality, (12.51) . (14.3p . (11.271) . and the character 
of curlv (by Lemma ETTl curlv G L^(M^) D L°°(M?)), we find that the last two 
terms in the right-hand side integral give o(l) terms. To bound the contribution 


< Cllv - v||lo 


\Vus — iUeN^v\ \ curlv| 
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of the first term, we split the integral over using the balls Bi given by 
Proposition 14.41 as follows : 


Vur — iUrNpYW curlvl 


< / |Vrte — iWeA^gvIl curl v| + / |Vtt£ — meA^evH curlv 

JuiBi 


< C\\Vue-iusN^y\\L 2 ,oo\UiBi\+i / \Vus - \ || curl v||£, 2 (r 2 ), 


/R2\UiR 


where we used ()4.12l) and the boundedness of curlv. In view of (I4.16p . (I4.18p . 
(j4.3p and item 1 of Proposition 14.41 we deduce that 

[ \Vu,- zu^A^evlI curlvl < Ce-^^N, + CN^ 

Jr^ 

and inserting into (j4.29l) and using item 2 of Lemma 14.71 we deduce that 


/ {je - iVeV) • (v - v) curl V 

4R2 


<o(l). 


For the second term of we apply item 3 of Lemma 14.51 with ^ = v — v to 
obtain 



27r^dj(v - v)(ai) +o(l)||v -vUc-o.t 
i 

< 27r ^ |dj|||v - v||l°o + o(l)||v - v||co.7 
i 


o(l), 


where we used item 2 of Lemma 14.71 and the fact that |(ii| < CN^ for these 
balls. We deduce that = o(l) and conclude that Id + Im + Iv = o{N^). 


4.3.3. The dominant terms. The term Is can easily be treated with the help of 
(j2.17p . Using that from Lemma 14.71 Vv vanishes outside of the balls Bi given 
by Proposition 14.41 and is bounded otherwise by a constant depending on v, and 
using (I2.17p . (I4.16P and Lemma 12.31 we may write for each t <T^, 



< 2||Vv||lcx> f \VUe - iUeNev\‘^ + ^{1 - \Us\‘^f + N^\v\‘^\l - \Us 

4R2\UiBi 

< C{£i;{ui;{t)) - vriVellogel) + o{N^). 


We thus conclude that 


|/s| < 



vrA^ellog 



+ o(iV2). 


(4.30) 
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For the term In we replace v*- by v’- + (v —v)-*-, and using Young’s inequality, 
we write 


(4.31) 
Id < 


N, 


|logel 


+ ||y - v||lo 


4 / / \dtUe-iUsN^4>f + 2 

^ Jo JR2 Jq J]r2 

N, 



\loge\ 



0 ./R2 


{VUe - iUeN^v) ■ V 
\dtUe - iUsN^c/)f + [ [ |VUe - iUeNeY\" 

Jo Jm? 


+ 


f f Ne\dtUe - iUeNs(l)\\divv - 4>\. 

Jo iR2 I log e I 


We claim that 


(4.32) 



Ne\dtUe 


iUsN^cpW div V 


I logs 


</>! 



iUeN^4>\‘^ + o(iV|). 


Indeed, either <C |loge|, in which case divv = 0 and </> = p and the result 
follows from the Cauchy-Schwarz inequality after inserting a factor 


1 

T^e 



1/4 


and the character of p; or llogej/Yg A and we may use (I4.28p and the 
character of divv to conclude the same. On the other hand, by (|4.3p and 
item 2 of Lemma 14.71 we have 


(4.33) ||v —v||loo 


E—^ ( f [ \dtUe — iUeN^Jil^ + [ [ IVUg — iUeN^vl 

|loge| \Jo U' ' ' Jo Jr^' ' 

< Ce~^ , [ [ \dtUe - iUeN^cpl^ + o{l). 
Jo JR 2 


|loge| 


We next distinguish two cases: 

Case 1: the case where 

[ [ \dtUe - iUeNi;4>f < 20 

Jo ./R 2 

By (14.3h this implies that 


|(Vtt£ - iUeN^v) ■ v|^. 


(4.34) f f \dtUe-iUeNs(j)\‘^<20\\Y\\L^{TTNe\loge\+N^). 

Jo 4r2 

It follows, together with (14.31) and (11.271) that < OiVellogel where is as 
in Proposition 14.81 From Lemma 12.11 we have that dtv is uniformly bounded 
while V G in space, hence v is Lipschitz in space-time, so we may apply 


















MEAN FIELD LIMITS OF GINZBURG-LANDAU 


37 


Proposition 14.81 with r = t, A = 2 and X = v to obtain 

N, 


(4.35) Iv < 


|loge| V 2 


- [ [ \dtUe 

2 Jo U' ^ ^ 




+ 2 


/V 

Jo JR 2 


Ue - iUsN^v) • vp + o{NJ). 


Case 2: the case where 


- iUeJ^eW >20 / |(VUe - iUeNeV) ■ v| 

Jo Jm? 


f j \dtUe 

Jo Jm? 

We may rewrite that condition as 

(4.36) i / f \dtUe - iUi;Ne(l)\‘^+ 4 [ f |(Vn£ - zUeiVev) • vp 
4 Jo ./R 2 Jo ./R 2 

[ [ \dtu^ - iUsNsCpl"^ + 2 [ [ KVue - meiVev) • vp. 

Jo Jk.^ Jo Jr^ 


, 1 1 

- * 4 10 


We note that in that situation, thanks to Lemma KT\ and the character 
of (j), we have < 2 fofjg ^2 \dtUe — iUsN^(j)\‘^ < CNj\logef, where is as in 

Proposition 14.81 Choosing Mi; = we may then apply that proposition 

with r = t, A = 4 and X = v, and combining the result with 

(j4.36D . we are led to 


N 1 1 C 

(4.37) + I \dtus - 


\loge\ 4 10 Jo 7 r 2 


+ / f |(V«,-m,iV,v)-v|2+o^ 


I log e I Jo 7 r 2 

This implies that 
(4.38) 


|loge| 


r/ \dtUe 

Jo ./R2 


- )+o(l) 


Iv < 




|loge| 


/ \dtUe - iUeNe(p\^ + 2 [ f | (V 

2 Jo Jr^ Jo Jr^ 

f f \dtUe 

Jo Jr^ 


Ue - iUeNeV) ■ v| 


1 Ne 
8 |loge 


- iUeNe(f)f + o(l). 


Returning to the general situation, we may now combine in the first case 
(j4.3ip . ()4.32l) . ()4.33l) . (I4.34p and (I4.35p . and in the second case (|4.31l) . (I4.32p . 
(j4.33p and (|4.38l) . Noticing an exact recombination of the terms 


1 /•£ /* pt 

i / / ,^ . ,|2 


\dtUe - iUeNe(p\ +2 

Jo Jr^ Jo iR2 

1 



(Vtte - iUeNeV) ■ V* 



+ o/ / \dtUe - iUeNeC/)J + 2 / \{VUe - iUeNeV) ■ v\ 

2 Jo JR 2 Jo JR2 
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we obtain in both cases that 

2N /"* f 

(4.39) Id + Iv<t. - ^ / |Vit£ - meiVeVpIvI^ + o(iV|). 

I log ^ I Jo Jr^ 

On the other hand, from (|4.19l) applied with ^ = |vp, we have 

2Ne I I ^ | |V-He - iUeNeV\‘^\v\‘^ - |vp/ie 

JR2 I logs I 

< (l|v||i- + “ ^^e|loge|) + o{N^), 

so using that N^- < 0(|loge|) and combining with (|4.39D . we obtain that 

(4.40) Iv + Id+Ie<C f (Seiueis)) - 7riV,|loge|) ds + o{N^). 

Jo 

Let us point out that this is the only place in the proof where we really are 
limited to the situation where < 0(|loge|). 

4.3.4. The Gronwall argument. Combining ()4.40p with ()4.30l) and the result on 
the negligible terms, we are led to 

£e{ue{t)) -Seiu^) <C ( (Te(ue(s)) - vriVellog e|) ds + o(N^) 

Jo 

and this holds for any t <T^. 

In view of the assumption on the initial data, Gronwall’s lemma immediately 
yields that 

(4.41) Se{ue{t),t) < viA^ellogel +o{N‘^) 

for all t < T^. Thus we must have Tg = min(l,T) and we may extend the 
argument up to time T to obtain that ()4.4ip holds until T. This proves the 
hrst assertion of Theorem \2\ 


4.3.5. The convergence result. To conclude the proof of Theorem [2l there re¬ 
mains to check that this implies that y in L^^^(M^) for p < 2. In 

view of (j4.16p . (I4.4ip implies that for every t < T, 

(4.42) [ \Vue- iu,N,v\^ < o{N^), 


hence for any ball Bji centered at the origin and any 1 < p < 2, by Holder’s 
inequality. 


(4.43) 


L 


BR\UiBi 


\VUe - iu^N^v\P 


< o{NP). 


Meanwhile, by (I4.18P and using the embedding of into L'^{Br) for 

any g < 2, by Holder’s inequality and item 1 of Proposition 14.41 we deduce that 
for any p < q < 2, 



iUeN^Y\P < o(l). 


Combining this with (j4.43ll . we conclude that ■^{Vu^—iUi.N^v) —)■ 0 in 
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In the case —)■ A, with (|2.5I) and the upper bound (j4.41l) . we have in 

addition that -^(Vue — iu^N^v) is bounded in It thus has a weak 

limit /, up to extraction of a subsequence £k- As in [SS51 p. 151], letting 


denote Uj-Bj for each given e, since r* < we may extract a further 

subsequence such that An '■= yJk>rS^£t, has Lebesgue measure tending to 0 as 
n —)• 00 . For any hxed n, by weak convergence we have 


lim inf 
k—^oo 




jVtte - iUsN^vf 

N! 


> lim inf 
k—^oo 


\VUe - iUsN^vf 

N! 


> 


/ 

Jr^\A 


l/l^ 


but the left-hand side is equal to 0 by (I4.42p . so letting n ^ 00 , we deduce 
that / must be 0. Since this is true for any subsequence, we conclude that 
■^(Vue — iu^N^v) converges weakly in to 0. 

The appropriate convergence of je/N^ is then deduced by (|2.6H (|2.7H and 
(j2.8|) , and this concludes the proof of Theorem [5J 


Remark 4.10. In order to treat the mixed flow or complex Ginzburg-Landau 
case, the computations are very similar and shown in Appendix [3 below in 
the gauge case. One should multiply the result of Lemma \2.4\ by 2fl\ + 
instead 0 /v*-, and use f = p or Adivv (respectively), and if = 13(f) — jvp. A 
supplementary error term inO{l3 fj ^2 14-(v—v)) appears, which can be controlled 
only by an estimate on /jg 2 IK] and leads to the extra condition ^ log jlogej. 


Appendix A. Proof of Proposition 14.81 

As already mentioned, the result is a quantitative version of the “product 
estimate” of |SS3] . It also needs to be adapted to the case of an infinite domain, 
which we do by a localization procedure based on a partition of unity. 

As in [SS3] we view things in three dimensions where the first dimension is 
time and the last two are spatial dimensions. By analogy with a gauge, we 
introduce the vector-field in three-space 

(A.l) Ae = Nfl(f>,v), 

whose hrst coordinate is N^<f> and whose last two coordinates are those of N^v. 
Equivalently, we can identify A^ with a 1-form. We also note that 

(A.2) ^curlAe = (curlv,(9tvi -di(j),dtV 2 -d 2 cf). 

We then define the 2-form 

(A.3) Je = d{{dUe - iUeA,;,iUe) + A,;) , 

where d corresponds to the differential in three-space. 

Lemma A.l. Identifying a spatial vector-field X with the 2-form X 2 dt Adxi -|- 
Xidt A dx 2 , we have 

(A.4) = fledxi A dx2 + Ve + {1- \Ue\'^)NfldtY - V(/>), 

where fi^ is as in ()2.1ip and as in ()2.12p . 
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Proof. By definition, = J^dxi A dx 2 + J^dt A dxi + J^dt A dx 2 , where for 
k = l,2, 

f J* = curl ((VWe - iUeNeV, iUe) + iVeV) = fl^ 

\ = dt{{dkUe - iUeNeVk, iuf) + NeVk) “ dk{{dtUe - iUeN^Cp, iUe) + NeCf) 

To obtain the expression of Jg, it thus suffices to compute 

dtiiVue - iUeN^v, iUe) + N^v) - V{{dtUe - iUsNef), iUe) + NeCf) 

= {VUe - iUsN^V, idtUe) + N^dtV + {dt{VUe - iUeNeV),iUe) - {VdtUciUe) 

- {dtUe,iVu,) +Vi{\u,\^ - 1)N,(P) 

= 2{VUe - iUeN^V, idtUs) + NsV{Ue,dtUs) + (1 - \Ue\‘^)N^dtV - NeY{dtUe, Ue) 
+ iV,V((|u,|2-l)0) 

= K - 2Ne<P{VUe,U,) + (1 - \Ue\^)NedtY + N,V{{\Ue\^ - !)(/.) 
where we used (I2.12jl . and this yields the result. □ 

We work in the space-time slab [0, r] x We consider X (here a spatial 
vector held, depending on time) and Y (here Y = et the unit vector of the time 
coordinate) two vector helds on [0,r] x In order to reduce ourselves to the 
situation where X is locally constant, we use a partition of unity at a small 
scale: let be as in (I4.2ip and let us consider a covering of [0, r] x by 
balls of radius 2Me centered at points of Mg and let {i^fcjfceN be an 

indexation of this sequence of balls and {xfcjfceN a partition of unity associated 
to this covering (which we observe has bounded overlap) such that X^fceN Xfc = 1 

and llVxfcll < . For each fe G N, let then Xk be the average of X in Dk- 

Then, working only in Dk, without loss of generality, we can assume that Xk 
is aligned with the hrst space coordinate vector ei, with ( 64 , 61 , 62 ) forming an 
orthonormal frame and the coordinates in that frame being denoted by (t, w, a). 
We will assume hrst that Xk 7 ^ 0. Let us dehne for each k, a the set 

= {(t,w)\{t,w,cr) G Dk}, 

which is a slice of Dk (hence a two-dimensional ball). Let us write for 

<^£( 61 , 64 ) restricted to ^k,a- In other words, by (IA.3p . if ^ is a smooth test- 
function on we have 

(A.5) / ^ A Js,k,<7 ~ / A {{du^ iu^A^, iu^) Y A^') 

where d denotes the differential in the slice 

We let Qk be the constant metric on 0,k^a dehned by fffc(ei, 6 i) = A/\Xk\, 
9kiet,et) = 1 /A and gk{ei,et) = 0 with A > 1 given. 

We then apply the ball construction method in each set ^k,cT- Instead of 
constructing balls for the hat metric, we construct geodesic balls for the metric 
associated to gk, i.e. here, ellipses. 
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Lemma A.2. Let C be as above and denote ^ = {x G Clk,a\dist {x, dQk,^) 

e}. Assume that 


Fe,k,a '■=\ j {\dtUe - iUeNe(j)\^ + \diUe - iUeN^Yi\^ + ^(1 - \Ue\^f 
^k,u 

+ A2(|V0|2 + |atv|2)) <M, 

with Mg as in (I4.2ip . Then if e is small enough, there exists a finite collection 
of disjoint closed balls {Bi} for the metric gk of centers a, and radii ri such 
that 
( 1 ) 

(2) iJiBi covers {||iig(x)| - 1| > ^} n 9.%^^. 

(3) Writing di = deg{u£,dBi) if Bi C ^ and di = 0 otherwise, we have 
for each i. 


(A. 6 ) ^ j^^(^A\Xkf\diU£ - iUeN^Yil'^ + ^\dtUe - iUsN^fif 

+ 2M-2A2(A|Afc|2|V<^|2 + ]-|5iv|2)) 

>7r\di\ (|loge| - ClogMe). 


(4) Letting ^i£^k,u = 27r djcJai; we have for any 0 < 7 < 1, and any 


/ 


C ^ Je^k^a ^f^£,k,c 


<cm\co,A^Mfi^yF£^k, 


Proof. The first 3 items are a rewriting of |SS3[ Proposition IV. 2], itself based 
on the ball construction, that needs to be adapted to the case of the nonstandard 
metric. As in [SS31 Proposition IV. 2] we start by noting, via the co-area 
formula, that there exists with ^M“^ < mg < such that setting tv := 
{|ug| < 1—mg} has perimeter for the gk metric bounded by CeM^. We may then 
apply |SS5[ Proposition 4.3] outside that set to V£{t,w) = ■i}}^(\/A|Xfc|r(;ei + 

and A£{t,w) = Ag(\/A|Afc|u;ei + -^Wt) restricted to the slice, with 

initial radius tq = CeM^ and final radius ri = . This yields a collection of 

disjoint closed balls Bi with sum of radii bounded by and such that 


1 

2 



+ Mg curl Ag|^ 


> TT^ |di|(|loge 

i 


Clog Mg). 


Making the change of variables x = '/h\Xk\w and s = we obtain balls Bi, 

the images of the BiS by the change of variable, which are geodesic balls for 
the metric gk and whose sum of the radii is bounded by (since A > 1 ); 

and inserting (IA.2p and using (14.211) we obtain (IA. 6 I) . We note that the fact 
that the domain size also depends on e does not create any problem in applying 
that proof. 

Item 4 is a consequence of Lemma 14.31 adapted to the present setting with 
differential forms, replacing Vug — MgAg by du^ — iu^A^ and using again (IA.2I) . 
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□ 


We now proceed as in the proof of [SS3j . We set Ve,k,cF to be the of 

Lemma lA.21 (item 4) if the assumption is verified, and 0 if not. We 

note that 


(A.7) 


\J. 


e,k,a t^e,k,a\\(cy{nk. 


< CAM- 


■1/2 P 


e,k,i 


is true in all cases. Indeed, either in which case the result is true 

by item 4 of Lemma IA.2I since Mg > M“^, or = 0 in which case, 

for any ^ G starting from (IA.5h and writing |(Vtie — iA^Ue,iUe)\ < 

iVug — iA^Usl + l|iiel ~ ~ we obtain with the Cauchy-Schwarz 

inequality, using the boundedness of ^lk,cT, 



+ l^el + |1 



< C'll V^llioo < Mg ^^^Fe,k,a, 

for £ small enough. But since F^^k ,(7 > we have y^F^^k,a + ^F^^k.a < 

2M£ and thus we find that (jA.7p holds as well. By (IA. 6 p . we also 

have that 


/ C 1 

|i^£,fc,(T| ^ lloge] 

Next, we choose rj a function depending only on time, vanishing at 0 and r, 
such that 77 = 1 in [Mg — Ms if that interval is not empty, and affine 
otherwise. By construction 


(A.9) WrjWco.i < CMl/\ llmllcoa < 


We may now write that 


(A.IO) 


In 


k,CT 


VXk 


’\diUe - MeiVeVip + -^\dtU^ - 
+ 2Mg-2iV|(A|Afcp|V</.p + ^|9ivp)) 

> (|loge| - Clog Me) / {rjxk - 

Indeed, if we are in a slice where ye^k.a = 0? this is trivially true. If not, we 
apply pA. 6 p and obtain 


\ / T^\{A\Xk\^\diUe-iUeNeYl\^ + \\dtUe-iUeNe(l)\^ 

^ |Afc| V A 

+ 2M-^Nl{A\Xk\^\Vct>\^ + t|5,v|2)) 

> 271^] |di|nun(? 7 Xfc)(|loge| -ClogMe)). 

I 

Inserting then in ing. ■ qxk > { vXk)i ai) - CArjllT^x^Uco,! and PA.9I) . and using 
item 1 of Lemma IA.21 yields (lA.lOp . 
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Next, we integrate (lA.lOp with respect to a, and combine with (1X71) and 
(lA.8h to obtain 


(A.ll) 



+ ^\dtUe 




+ 2M-2iV2(A|Xfc|2|V0|2 + ^|5iV|2)) 

> (|loge| - ClogMe) / VXkJeA^ 

7]R2 X [0,r] 

- C’max(-^, A|Afc|)A2Mi/4-iF,,fc + CM^/^\loge\AM-^/^ 


with = := f F^^k,ad(T. Moreover, by (jA.4p . (since we assumed is along 

the direction ei) we have 




e^k,<j 






K • Xfc + X,(l - \Ue\^){dtV - VcP) ■ Xk 


so we may bound 


/ 

JDk 


r]XkXeil - \Ue\‘^)idtV - V(t)) ■ Xk 

< C||X||ic.iV,||l - Kl^WmnJdtv - 


< C||X||iooAr^eF,,fc. 


Inserting into (lA.llh and multiplying by |Xfc|, we may write 


VXk(^A\Xk\‘^\dxUe - m^X^vp + ^ \dtUe - iu^Ne4>\‘^ 

+ 2M-2iv2(A|Xfc|2|V</.p + ^|9iv|2)) 

> (|loge| - (7 log Me) f rjXkVe ■ Xk 
JDk 

c(i,A|Xfc|2)A2M-3/4 + (AM-V4|ioge| +eiV,)||x||iocFe,fc, 


— max( 


and we note that this holds as well if Xk = 0. We may next replace Xk by 
X in the left-hand side and the J" 14 ' term, and using that |X — X^l < 
CM4^/^||X|| qo,i in Dk, the error thus created is bounded above by 


X4-^/^A|loge|||X|4oc(l + \\X\\co,i)F, 


e^k 
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where we have used that by definition of V^, we have |^| < We may 
thus absorb this error into the others and write (since A > 1 ) 


/ VXk 
JDk 


’\dxUe - itteiVeVp + ^ \dtUe - iUeNs(f)\^ ^ 


+ 2M-2iV2(A|x|2|V0|2 + l|aiv|2)) 

> (|log e| - C log A/e) [ r^XkVe ■ X 

JDk 

- C\\X\\l^ (a 3(1 + ||X|bo.i)M,-V8 + 


Summing over k, using that YlkXk = 1 in [0 ,t] x the finite overlap of the 
covering, the fact that dtv and Vcj) S A^([0, r], by Lemma [XT] and (12.3p . 

we are led to 


'0 /R2 


r] ( A\X\‘^\dj;Ue - meiVeVp + F\dtUe - me(/)| 


A 


> (|log e\-C log Me) [ f T]Ve-X 

Jo /r 2 

-C\\X\\l^ (A3(l + ||X||co.i)M-i/«+eiVe) {Fe + N^) ■ 

Moreover, by choice of i] and definition of V^, we have 

r f {l-v)X-Ve 
Jo /R 2 

< 2||X|koo / / \dtUe — iUeNe(j)\\VUe — iUeNeVl 

< C\\X\\loo IFe( sup £eit) + o{l))M-F\ 

V ^ie[0,r] ^ 

where we have used the Cauchy-Schwarz inequality and (12.51) . Combining the 
last two relations, we deduce the desired result. 


Appendix B. Existence and uniqueness for (I1.14p 


We recall that the definitions of Holder spaces that we use are at the be¬ 
ginning of Section [2j For the sake of generality we study (|1.14l) with arbitrary 
a > 0 and /3 > 0 such that = 1, and A > 0. We denote by a; = curlv 

and d = divv. We note that if v solves ()1.13l) then (w,/) solve 


(B.l) 



2 div (fJvu! + aw^uj) 

—Ad -|- 2 div — avoj). 

a 


Theorem 3. Assume A > 0, a > 0 and 0 ^ + 13“^ = 1. Assume v(0) is such 
that ^a;(0) is a probability measure which also belongs to and d(0) G 

C'^ n for some 0 < 7 < 1 and some 1 < p < 2. Then there exists a 

unique local in time solution to (I1.13P on some interval [0,T], T > 0, which is 
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such that V £ for any < 7 . Moreover, we have v{t)—v{0) G 

L°°([ 0 ,r],L 2 (M 2 ))^ d^y g Lmo,T],L^{R^)), d £ L°^{[0,T], LPiR?)) and ^uj{t) 
is a probability measure for every t £ [0,T]. 

We start with some preliminary results. 

Lemma B.l. Let v be a vector field in L°°([0, T], cjq £ C'^, 0 < 7 < 

1, / G C^) then the equation 

/R ON / dtuj = div {vu) + / 

^ ^ ^ 1 a;( 0 ) = cao 

has a unique solution, and it holds that for some (7 > 0, Cq > 1, 

pt \ pt 


(B.3) ||w(t)||ii(R2) < ||wo||i,i(r 2)+exp ||Vu(s)||loc dsj ||/(s)||l°° ds 

and for any — 1 < u < 7 , 

(B.4) 

||w(t)||c'" < Co ^llwollc'" + I ||/(s)||c<T ds + ||a;o||i,°° J || divu(s)||c<T ds 

X exp (c ||Vu(s)|| ds ) . 

Proof. One may rewrite the equation as 

dtUJ = u • Vw + w div V + f. 

Then, by propagation along characteristics, we obtain first 




1,1 + 


I II/WII exp J II div u(s) 11^00 ds, 


second. 


|w(i)||L°° < ^||wo||l°° + J \\f{s)\\L-^ ds 


exp / II div u(s) 11^00 ds, 

Jo 


and third (IB.4p for <7 = 7 follows by a Gronwall argument. For general u < 7 , 
one can proceed as in |BCDl Chap. 3]. 

□ 


The next lemma about the regularizing effect of the heat equation can be 
found in |Ch2[ Proposition 2.1] (applied with p = 00 , p = 00 and noting that 
is the same as C^). 

Lemma B.2. If g £ T°°([0, T], and uq £ then the 

equation 


(B.5) 


J dfU = vAu + g 

\ u(0) = Uq 


has a unique solution which is in L°°([0, T], C''>'nL^(M^)) and T^([0, T], H^{E?)) 
and for any <7 < 7 , 


(B.6) 


\\u 


li:,°°([o,r],C‘^) < Co ^Ikollc*^ + 'y^lls'llL°°([o,r],c-i+'^)^ • 
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We note that the fact that u G L'^{[0,T], (B?)) comes from the fact that 

g G L^([0,T],and the regularizing effect of the heat equation. 


Lemma B.3. Assume u solves in [0,T] 
(B.7) 


dfU = Au + div / 
u{0) = Uq 


then, if q < p and ^ ~ ^ 5 > 0, we have for any t G [0, T], 


< 


LP + Cp^qfP 9 2 


L°°([0,T],L9(R2))- 


Proof. We follow |Du[ Lemma 2.3]. Using Duhamel’s formula, we may write 

u{t, •) = G{t, •) * no + w{t, •) 

with 

w{t,-) = J j Tt-s{x - y)f{s,y)dsdy 

_IXI ^ 2 

where G{t,x) = and rt(x) = -dxG{t,x) = By Young’s 


inequality for convolutions, we have 


lk(^)llLP(M2) < ||wo||lp(R2)I|G(L + II'“^(^)IIlp(R2) 

< lko||LP(R2) + ||^i’(i)||LP(R2)- 

We turn to studying w. We may write with Holder’s inequality, 
t / \ ^ / 

Tt-s{,x-y)i\f{s,y)fidy ) ‘ ds 


\w{t,x)\< j (^J Tt-six - y^^ dy'^ 


with q' such that 1/q + 1/q' = 1, and hence with Holder’s inequality again, if 
Q<P, 


I'W^(^)IIlp(R2) 

<r 


■ 9_ \ 9 ' 


Tt-s{x - y) 2 \f{s,y)\'^ dy] dx ] ds 


< 


l^ijrl.dyy 


(/ Tt-six - y)^yf{s,y)y dy){ I |/(s, y)]''dy)''? ^dx]' ds 


< 


Q_ \ q' 


|/(s,y)|'?dy 


q p 


P 


t—s 


|/(s,y)|'?dy ds 


10 




\f{s,y)ydy 


P \ P 

r/_J ds 


where for the passage from the second to third line we used Young’s inequality 
for convolutions. We may thus write 

11/2 

11-^5 1 

'0 


|R^(t)||LP(R2) < ||/||L-([0,r],M(R2)) / \\^s\\]^,^^2^^2)\\Ps\\lp/2(^2)ds. 
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But computing explicitly, we find ||rs||^T-(]R 2 ) = CrS 3 / 2 +i/r deduce 

that 

II^(^)IIlp(R 2) < C'p,gll/llL°°([0,T],L9(R2))t2 P 9, 

and the result follows. □ 


Finally, we will need the following potential estimates: 


Lemma B.4. Let u e C'^ n LP{E?) with 0 < 7 < 1 and 1 < p < 2. Then 
VA“^u, where is meant as the convolution with —^log, is well defined 
and 


(B.8) ^^IIc1'7(R2) < C'l(lkllc^(R 2 ) + II^IIlp(R 2 ))- 

Proof. Setting v = VA“^u we may write 

and one may check that this integral is convergent with 


(B.9) 

|u(x)| < (711^11^00(^2) 


'|y—a:|<l 1^ 2/1 


dy + \\u\\lp 


1/p' 


\y-x\>i \x-y\p' J 

^ C'(II^IIl°°(r2) + II^^IIlp 


with 1/p + 1/p^ = 1. Let then w be such that v = Vm, and thus Aw = u. For 
any x € M^, by Schauder estimates for elliptic equations and (IB.91) . we have 

||V'«;||c'i,7(s(x,i)) < C'(l|Vu'||i:,oo(B(2._2)) + \\ u \\ c '' i { b { x , 2 ))) 

^ C'(||u||j;^p(r 2) + ||tf||c'7(R2))- 

Taking the sup over x yields the result. □ 


We now turn to the proof of Theorem [3] for which we set up an iterative 
scheme. 

Let Vo = v(0) andcjo = curl vq, do = divvQ. Given v^ andcun = curlv^jd^ = 
divvn, we want to solve 


(B.IO) 


Let 

where 


dtUJn+i = 2div ((/3vn + av^)wn+i) 

, dtdn+l = ^Adn +1 + 2 div ((/3v/; - aVn)Wn) 
a;„+i(0) = a;(0) 

_ d„+i(0) = d(0). 

tn = sup{t,i(:„(t) < 2CiCqKq} 


Kn{t) = ||a;„(t)||c7 + \\uJn{t)\\Li + \\dn{t)\\c^ + ||dn(t)||LP + ||Vn(t) HcLt 

Kq = ||a;(0)||ii(R2) + ||a;(0)||c7 + ||d(0)||LP + ||d(0)||c7, 

and Co is the maximum of the constants in ()B.4D . ()B. 6 D . and Ci the maximum 
of the constant in (|B.8|1 and 1. 
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Let US show by induction that (IB.lOp is solvable and there exists Tq > 0 
independent of n such that tn > Tq for all n. 

For n = 0, this statement is true by assumption. Assume it is true for re. 
Then in view of Lemma IB.11 we have that the equation (IB.lOp for re + 1 has a 
solution Ldn+i G L°°{[0,T],C'^ n L^), with, if T < Tq < tn, 

(B.ll) 

\\^n+l\\L°°{lO,T],C-y) 

< C'o(^||w(0)||c7 + ||a;(0)||Loo J \\vn\\c^.jis) ds^ exp j ||v„||ci (s) 

< Co (||u;(0)||c7 + \\u;{0)\\L^2CiCoTKo)exp{2CiCQCTKo ). 

Also it is straightforward to see by integrating the equation and since con+i G 
by (IB.3P that ^uin+i remains a probability measure. 

Similarly, Lemma IB. 2 1 yields the existence of a solution dn+i G L°°([0, T],C'^) 
with, if T <Tq <tn, 

(B.12) ||d„+i||^oo([o,T] ,CT') < co{\\dm\c^+ cVt\\ YnWrillcr) 

<Co{\\dm\c.+cVfCiClKi), 

where C depends only on a and A; and by Lemma [B]3] applied with q = 1, 

1—1 

(B.13) ||d„+i ll^^ooQo,T],LP) < \\d{0)\\LP + CTp ^ ||Vn<Vn+l||L°o([0,T],I,l) 

< \\dm\Lp+CT'p-^CoCiKo. 

We then let 

(B.14) Vn+l = VA~^dn+i - V-^A'^Wn+i- 

By Lemma IB.41 this is well-dehned and 

(B.15) ||Vn+i 11 ( 71,7 < Cl (||dn+i||(77 + ||dn+l||LP + ||^Vn+l ||c")' + ll‘^n+l||Ll) • 

In view of (jB.lip , (IB.12P , (IB.13P , and (IB.lSp we then deduce that if Tq is chosen 
small enough (depending on Kq and the various constants), then we will find 
that tn+i > Tq. The desired result is thus proved by induction. 

Let us now show that {w„} and {dn} are Cauchy sequences in C~^~^'^. By 
subtracting the equations for re and re + 1 we have 
(B.16) 

dt{uJn+l - Wn) = 2div (^{(3Vn + a\^){uJn+l - Wn) 

+(/3(vn - v„_i) + a(v„ - v„_i)-^)a;„^ 

dt{,dn+i dn) = dn) 

+2div ((/3v;|; - av„)(a;„ - ujn-i) + {/3{vn - Vn-i)"^ - a{vn - v„_i))a;„) . 
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Since = a;(0) for all n, applying (IB.4h with cr = —1 + 7 , and the result of 

the previous step, we have 

||Wn+l — l^n||L°°([0,r],C‘^) 

< e^fo IIVv„|Uoo ds f II ((/3(v„ - Vn-i) + a{vn - v„_i)-^a;n)||c- ds 

Jo 

— Ce^'^T\\{Vn — Vn-i)uJri\\L^{[0,T],C^+’^) 

< Ce^^’^T (||a;„ — 0Jn-l\\L°°{[0,T],C'’) + \\dn — '^n-l||L°o([0,r],C‘^)) • 

Similarly, in view of (IB. 6 p . 

\\dn+l — C^n||L°°([0,T],C'^) < CVT (^\\VnioJn “||L°°([0,r],C‘^) 

+ II (Vn — Vn-l)Wn||L°°([0,T],C'^)) 

< CVt (||a;„ — w„_i||x,oo([o,r],C‘^) + \\dn — i^n-i||L°°([o,T],C'^)) > 

if (T < — 1 + 7 . We deduce by standard arguments that {oJn} and {dn} form 
a Cauchy sequence in L°°([0, T], C"^) if T is taken small enough. By inter¬ 
polation {vn} is a Cauchy sequence in , for 7 ' < 7 (resp. {oon} in C'^ 
and {dn\ in C'^ ). The limits v, oj^d will obviously solve (jB.lh . oj will be in 
L°° ( [0 , T], (M2)), d in L°° ( [0, T], LP (M2)) n ( [0 , T], (M2)) by Lemmas EJ 

and lB.3[ and (from (IB.14M 

(B.17) v = VA-^d-V^A-^w, 

with V G L°°([0, T], Taking the time derivative of ()B.17p and using (jB.ip . 

it is standard to deduce that v must solve (I1.14p . By integrating the equation, 
using that v is bounded and oj G we also have that remains a 

probability density. 


Next, we prove that v —v(0) remains in Using (11.141) and integration 

by parts, for Q{x) = with 0 < 7 < 1 we compute 

- 7 - f C|v(t)— v(0)p = 2 [ C(v—v(0))-( — V div V-|-2/3v*-curlv — 2Q:vcurT 
dt Jm. 2 J^2 \a 

j C(divv)2 + ^ [ C(divv)(divv(0)) - ^ [ VC • (v - v(0)) div v 


a 


O JR2 


« dR2 


4q; f Clvl^w-|-Cv • v(0)a; — 4/3 f C^"*" • v(0)a; 
yR2 dR2 


and using Young’s inequality, the fact that |VC| < |CI that divv G we 
find 

d_ 
dt 


C|v-v(0)|2<c(/ C|divv(0)p+/ C|v(0)|2u;+/ C|v-v(0)|‘ 

VdR2 dR2 dR2 


with a constant C depending on a, /3 and A. Using Gronwall’s lemma, since 
Jjg 2 |w| is uniformly bounded and v( 0 ) is bounded, we deduce that fj ^2 C|v(t) — 
v(0)p < Ct, with Ct independent from 7 . Letting then 7 —)• 0, the desired claim 
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follows. 


To prove uniqueness, it suffices to apply the same idea of weak-strong unique¬ 
ness as used in the proof of the main results: if vi and V 2 are two solutions of 
(jl.ldp . we introduce the “relative” energy E{t) = ^ fj ^2 |Yi(t) — V 2 (t)P, and the 
relative stress-energy tensor 

T = (vi - V2) (8) (vi - V2) - ^|vi - V2|^/. 

By direct computation we have 


(B.18) 

—E= / (vi — V 2 ) (—V( divvi — divv 2 )-|-2/3(v]’‘curl vi — v^ curl V 2 ) 
at Jt ^2 V a 

— 2a (vi curlvi — V 2 curlv 2 )^ 

= [ - —I div (vi - V 2 )p-I-2/3(vi - V 2 ) • v^ curl (vi - V 2 ) 

JR2 Oi 

+ / —2a|vi — V 2 p curl vi — 2a(vi — V 2 ) • V 2 curl (vi — V 2 ). 
Jr'2 

and 

(B.19) divT = (vi — V 2 ) div vi -|- (vi — V 2 )''' curl (vi — V 2 ). 

Multiplying (IB.19I) by 2/3v2 -|- 2av^ and rearranging terms, we find 


/ (2/3v 2 -k 2av^) • divT = / (2/3v2 -k 2av^) • (vi - V 2 ) div (vi - V 2 ) 

Jr^ Jr^ 

-k / ( 2 av 2 - 2/3v^) • (vi - V 2 ) curl (vi - V 2 ), 
Jr'2 

and inserting into IjB.lSp and using one integration by parts, we obtain 


^E < - f — I div (vi - V 2 )p-k [ T : V(2/3v 2-k 2av^) 
hr Jr2 a J^2 

-k / (2/3 v 2-k 2av^) • (vi - V 2 ) div (vi - V 2 ). 
JR'^ 

We may next use the boundedness of V 2 and Vv 2 , bound |T| by 2E, and use the 
Cauchy-Schwarz inequality to absorb the last term into the first negative term 
plus a constant times E, to finally obtain a differential inequality of the form 
■^E < CE. The integrations by parts can easily be justified by using cut-off 
functions like XR rmd taking the limit, using the fact that Vj — v(0) G L‘^(M?) 
and div (vi — V 2 ) G We then conclude to uniqueness by Gronwall’s 

lemma. 


Appendix C. The gauge case 

In this appendix, we present in a formal manner the quantities that should 
be introduced and the computations which need to be followed to obtain the 
limiting dynamical laws in the two-dimensional case with gauge mentioned in 
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Section 11.3.41 of the introduction. We use the notation of that section, in par¬ 
ticular the starting point is (11.311) . In the rest of this appendix, we will not 
write down negligible terms, such as terms involving 1 — |ue|^. We will however 
present the computations in the general case of the mixed flow (ll.3ip . which 
can thus be used as a model for studying the mixed flow in the case without 
gauge. 

We start by introducing some notation: we set 

d<s> ■■= dt - i^e 


and 


(C.l) <^:=| 

We will also denote = 
(dropping the e) and 
Lemma 2.12 in m- 


p in cases leading to (jl.32l) 

^ div V in cases leading to (I1.33P . 

+ NeCp and dq>f = dt- A' = Ae + iVeV 

= curlAg. We define the velocity as in (2.43) and 


(C.2) 14 := + dt{iUe,VAeUe) - 

= 2{id,s,^Ue,VA^Ue) + (l^tel^ “ '^)Ee 

thus 

(C.3) dtje = Ee. 

The modulated velocity is then defined as 

(C.4) Ve = 2{idis>'Ue,VA'Ue) + (l^ep - l)Ee. 

The modulated energy is defined as 
(C.5) 

£6iUs,A) = l [ \VAeUe-iUeNeY\^ + \cUTlAe-N^h\^+ ^'^ ^ +{l-\Ue\^)lp 

2 Jm 2 2e 

with 

(C.6) = jdcj) — |vp. 

We will also use the fact that for Ug solution of (11.101) we have the relation 

(C.7) divje = NsHy, -r + iP)d<s>Us,me) ■ 

|loge| 

The stress energy tensor is now 


(C.8) {Se{u,A))kr.= 

{{dku - iAkU, diu - iAiu) - ^ “ I 

A direct computation shows that if u is sufficiently regular, we 

divS'£(tt,A) := 'Y^di{Se{u, A))u 

i 


curl Ap 
have 


^ki- 


= (Vau, - |rip)) 


curl A(V''‘ curl A + {iu, Vau))'^ 
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SO if Ue solves (11.1011 . we have 

Ol I 

(C.9) divS'£;(n£, Ae) = Ne{{j -r + VAe^e) + curl 

|loge| 

For simplicity, we will denote it by S^- The modulated stress tensor is now 
defined by 

(C.IO) 

{Se{u, A))ki = {dkU - iuAk - iuN^Vk, diu - iuAiu - mA^eV^) + (1 - |up)vfcVz 

^|Vau - mA^evI^ + (1 - |vl^ + ^(1 - \u\‘^f - I curl^ - iVehp^ 6ku 

and by direct computation 

(C.ll) div Se = div Se + N^v div v + curl v - div v 

- curl V - A^evdivje - curlje + V(ilV|h^ - N^hh^;). 


Combining (jC.Op and (1C.lip and we obtain 
(C.12) 

div,Se = N^{{ ^ + i(3)d^Ue, + aheE^ + A^^vdivv + curlv 

I ^ 1 

- div V - Nej^ curlv - Aievdivjg - AleV-^ curlje + - A^eh/le). 

Writing d<^Ue = d^'U^ + iu^N^cp a-^id VA^e = A'Ue + iUeN^v, and inserting 
(jC.2p and (1C.71) into (|C.12p . we are led to 


(C.13) divSe = |]^gg| {d<s>'UeyA'Ue) + crheE^ + V(^A^|h^ - N^hh^) 

+ ^NeVe - + A^^ curlv - A^^j^ curlv - NeV-^ curl je + o(l). 

Multiplying relation (|C.13p by 2/3v + 2av-*- yields 
[ (2/3v + 2av*-) • div5e 



+ 



+ 



+ ^NeVe - ■ (2/3v + 2av-^) 

2aA(|lvp curlv + N^j^- ■ (2/3v'‘ — 2q:v) curlv — 2aN^\v\‘^ curlj^ 
V(^A^eh^ — A^gh/ie) • (2;0v + 2av*-) — cr/igFig(2/3v*- — 2av) + o(l). 
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This allows to express fj ^2 ■ v so that splitting 14 • v as • v + /3^14 • v 

(since + /3^ = 1) and inserting the previous relation we may obtain 

(C.14) 

- [ NeVe-v= [ N^Ve ■ {-a'^v + a/3v-^) - N^f3‘^\v\‘^dt\Uef 

JR2 yiR2 

+ f {2Py+ 2av-^) ■ (-div Se{Ue) + A'Ue)] 

V |loge| J 

+ / 2aAif jvp curlv + ■ (2/3v*- — 2av) curl v — 2aNs\v\‘^ curing 

iR2 

+ [ V{-N^h.‘^ — N^hciiilA) ■ {2f3v + 2av^) — ahsE^(2Pv^ — 2av) + o{l). 
Jr2 2 

The analogue of Lemma 12.51 combined with (jl..S2h or (II.Tip is 
d f Of /* 

—£e{ue) = - + (t\E^\‘^ + / Ne{iUs, d^Ue) divv 

at Jr 2 I log e I Jr 2 

+ [ Nlv • dtY - Neje • dtv + N,{-V, - E,) ■ Y + N,dMN,h - h,) - N,WtK 

2r2 

+ [ l^edt ((1 - \Ue\'^){'>P - |vp)) . 

JR2 2 

Using ()1.32p or (ll.3.Sp . we may write 
[ N^y ■ dtv - NJ, ■ dtv 

JR2 

= J iVePA^gV — je) • + (—2av + 2/3v*-)(curlv + h) + V(/)^ . 

Next, we insert again (9$Ug = + iUsN^cp and write that 

(C.15) iVgV - jg = -V^(A^gh - hg) - a(iVgE - Eg) 

which holds by subtracting the equations (11.311) and (ll.32p - (ll.33p . and we insert 
(jC.14p and use an integration by parts to obtain 


d 

dt 


£e = f Sg : V(2/3v + 2 q : v '“) + iVgl4 • (— a^v + a/Sv*-) — N'g/S^lvpcltlng 

JR2 

- / 2aiVg|v|^ curl jg + N^cpdiY (iVgV - jg) 

4r2 

“ [ \ {d,^'U,VA'u) ■ (2/3v + 2av‘~) 

JR 2 |loge| |loge| 

+ [ iVg(-V^(iVgh - hg) - a{N,E - Eg)) • ((-2av + 2/3v^)h + E) 

./r2 

+ [ -a|Eg|2 - AT^Eg • V - ATg curl E(iVgh - /ig) + A^gh curl Eg 
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+ 


JR2 2 ! 


N^hhs) ■ (2/3v + 2aY^) — aheEe{2l5v^ — 2av) 


+ 





Next, we replace I 4 by I 4 + N^vdt\ue\‘^ — N^(j)V\u^\'^ obtained by comparing 
(|C.2|) and (jC.4|) . and use ()C.7h to reexpress divje;. By choice ()C.6p . the terms 
in factor of (9t(l —|rt£p) cancel, and the terms in factor of (j), div v and {d^Ue,iUe) 
formally cancel (up to small errors) as in the parabolic case. Then, we insert 

ah^Ei; ■ (2/3v'‘ — 2q;v) 

= ahe{E^ — N^E) ■ (2/3v'‘ — 2av) + Neh£{—v — V‘''h) • (2/3v'‘ — 2av) 

= ah£{E£ — N^E) ■ (2/3v’‘ — 2av) + 2aiVe|vp/i£ — N^heV^h. ■ (2/3v'“ — 2av). 

The term 2aN£\Y\^h£ gets grouped with 2aA^£|vp curl to form —2a|vp/i£. 
Combining all these elements, there remains 



Keeping only the last three lines, after some simplifications using again (|1.32[1 . 
and some integration by parts, these three lines can be rewritten as 


- he) 
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We thus finally obtain 
(C.16) 


d 

dt 

I 


Se= : V(2/3v + 2av'‘) + / • (—— / 2aN^\v\‘^fie 


Nr a 


\dq,’U, 


£|^ + Na'Us) ■ (2/3v + 2av-^) 


|loge| |loge| 

[ a{NrE-Er)i2av-2/3v^){Nrh-hr)+ [ -a\Er - NrE\^ + o{l). 

Jr2 Jr'2 

For the terms involving E we write with Young’s inequality 

[ a{NrE-Er){2aV-2/3Y^){Nrh-hr)<a [ \Er-NrE\^ + c[ |iV,h-h,|2 

Jr2 J^2 

and we deduce that the last line in (jC.lhp is bounded above by C J]g2 l-Yeh — 
hep, which we claim it itself bounded by a constant time the energy excess 
£r — 7rYe|loge|. Indeed, in writing down the analogue of Proposition 14.41 and 
14.161 one may replace f ^2 \ curlAg — iV^hp in (jC.SD by half of itself and still 
obtain the same optimal lower bounds for the energy in the balls, so we deduce 
that jg2 I curing — iV^hp must be bounded by the order of the energy excess. 
To finish, we use Young’s inequality again to bound 


/ 

JR- 


Nra 

|loge| 


{d<i,'UriNA'Ue) ■ (2/3v + 2av' 


< 


Nra 

|loge| 


19$/tie 


+ 2 


Wa'Uc ■ (/3v + av-^)p 


and the product estimate to formally bound 

[ NrYr ■ (-a^v + a/3v-^) 

JR2 


< 


Nea 

|loge| 


|5$/U£p + 2 


Wa'Us ■ (-av + /3v-^)p 


and using that av — /3v-*- = (/3v + av*-)-*- and |/3v + av-*-p = |vp, we see that 
these two relations add up to a left-hand side bounded by 


Nea 

|loge| 


\d^iUe\^ + 2 


[ \Na 

JR2 


'Ur\ Y 


which will recombine with — J ^2 [ Mge ] l^$'Wp+2Q;Yg|vp/rg into a term bounded 
by C{£e — vrY'ellog e|). Inserting into (IC.16h and replacing the use of v by 
that of V in the case with dissipation, we may then obtain a Gronwall relation 
■^£e < C{£e — vrY'ellogel) + o{N^) as in the case without gauge, and conclude 
in the same way. 
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